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ABSTRACT 

We obtain constraints on cosmological parameters from the spherically averaged 
redshift-space correlation function of the CMASS Data Release 9 (DR9) sample of 
the Baryonic Oscillation Spectroscopic Survey (BOSS). We combine this informa- 
tion with additional data from recent CMB, SN and BAO measurements. Our results 
show no significant evidence of deviations from the standard flat-ACDM model, whose 
basic parameters can be specified by fi,n = 0.285 ± 0.009, 100 f^b = 4.59 ± 0.09, 
ris = 0.961 ± 0.009, H^ = 69.4 ± O.Skms"! Mpc"^ and cts = 0.80 ± 0.02. The 
CMB+CMASS combination sets tight constraints on the curvature of the Universe, 
with J7fc — —0.0043 ± 0.0049, and the tensor-to-scalar amplitude ratio, for which we 
find r < 0.16 at the 95 per cent confidence level (CL). These data show a clear sig- 
nature of a deviation from scale-invariance also in the presence of tensor modes, with 
TT-s < 1 at the 99.7 per cent CL. We derive constraints on the fraction of massive 
neutrinos of f^, < 0.049 (95 per cent CL), implying a limit of ^ m^ < 0.51eV. We find 
no signature of a deviation from a cosmological constant from the combination of all 
datasets, with a constraint of wde = —1.033 ± 0.073 when this parameter is assumed 
time-independent, and no evidence of a departure from this value when it is allowed 
to evolve as wbe{cl) — wq + Wa{l — a). The achieved accuracy on our cosmological 
constraints is a clear demonstration of the constraining power of current cosmological 
observations. 

Key words: cosmological parameters, large scale structure of the universe 



1 INTRODUCTION 

In recent years, a wealth of precise cosmological obser- 
vations have been used to place tight constraints on 
the values of the fundamental cosmological parameters 
(e.g. iRiess et al." 1998; ' Perlmutter et a l'lOgQ*: 'Spcr gel et aT 



2003 



Riess et al.., 2004: iTegmark et al. 2004: S anchez et al 
Spergel et al . 2007; Ricss ct al. 2009; K omatsu et al 
2009; Komatsu 2010; Pcrcival et al 
20101: iReidet al.ll2010a: Ricss ct al. 2011; Blake ct al. 2011 



2006 



Montesano et all |2012| ). The unexpected conclusion from 
these studies is that we seem to live in a more complex 
and richer Universe than originally suspected; one which is 
currently undergoing a phase of accelerating expansion. Un- 
derstanding the origin of cosmic acceleration is one of the 
most outstanding problems in physics as it may hold the key 
to a true revolution in our understanding of the Universe. 

Within the context of general relativity, cosmic acceler- 
ation implies that the energy-density budget of the Universe 
is dominated by a dark energy component, which counter- 
acts the attractive force of gravity. A key parameter that 
can be used to characterize this component is the dark en- 
ergy equation of state woe, defined as the ratio of its pres- 
sure to density. In the standard ACDM model, dark energy 
can be described by a fixed equation of state specified by 
Woe = — 1, which can be interpreted as the quantum energy 
of the vacuum. However, a large variety of alternative models 
have been proposed, which predict different values of kjde 
and its possible evolu t ion with time (for a review see e.g. 
Peebles fc Ratrall2003l : iFrieman et ai]|2008l : iGott fc SlepianI 
20lj ). 

Measurements of the large-scale structure (LSS) of the 
Universe are expected to play a major role at shedding light 
on the causes of cosmic acceleration. The shape of the galaxy 
power spectrum, P{k), and its Fourier transform, the two- 
point correlation function £,(r), encode useful information 
which can be used to obtain robust constraints, not only 



on dark energy, but also on other important physical pa- 
rameters like neutrino masses, the cu rvature of the Uni- 
verse or details of inflationary physics (IPercival et al.l 2002 
Tegmarket"ai] l2004': Cole et al. 2005; Sanchez et al. 2006 

2009; Komatsu 20 1" 
Keisler et al.l 12011 



Spergel et all 



Percival et al 



2007; 
'2OIO: 



K omatsu et al 
R eid et al. | 



2010a : 



Blake et al.ll2011i : iMontesano et al.ll2012l 'l. A special feature 



of large-scale clustering measurements provides a powerful 
method to probe the expansion history of the Universe: the 
baryon acoustic oscillations (BAO). These are a series of 
small amplitude oscillations imprinted on the pow er spec- 
trum JEisenstein fc Hul [l99i : iMeiksin et "all 1 19991 ). which 
are analogous to the acoustic oscillations present in the 
cosmic microwave background (CMB) power spectrum. In 
the correlation function these are transformed into a sin- 
gle peak whose posit ion is related to the sound horizon 
at the drag redshift (jMatsubaral [200J) . As this scale can 
be calibrated to high precision from CMB observations, 
BAO measurements at different redshifts can be used as 
a standard ruler to measu r e the d i stanc e-redshift relation 
(|Blake fc Glazebrookl l2003l : iLindej |2003| ). The BAO fea- 
ture was first detected in the clustering pattern of the 
luminous red galax y (LRG) sample of t he Sloan Digital 
Sky S urvey (SDSS, lYork et al.l I2OO0I ) by lEisenstein etal] 



an d the Two-degree Field G alax y Redshift survey 
(2dFGRS, IColless et alli200ll . 120031 ') by Icole et all (|2005l ') 



and has been subsequently observed using a variety of 
datasets and techniques (jHiitsi 2006,: IPadmanabhan et al.l 
20071: iPercival et all 12007). boid: ICabre fc Gaztaiiagal l2009l: 



Gaztanaga et all 20091 : iKazin et a l.' '2010!: 'Beut ler et all 
2011; Bla ke et al.ll201ll : IHo et"aL| [2012: Sco ct al. 20l3). 

Driven by the potential of LSS observations for shed- 
ding light on the problem of the nature of dark energy, 
several ground-breaking galaxy surveys are currently being 
constructed or designed which will be substantially larger 
than their predecessors. The ongoing Baryonic Oscillation 
Spectroscopic Survey (BOSS. ISchlegel. White fc EisensteinI 
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[2009") is an example of these new surveys. BOSS is a part 
of SDSS-III (Eisc nstein et al. 20111 ) aimed at obtaining red- 
shifts for 1.5 X 10'' massive galaxies out to z = 0.7 over an 
area of 10,000deg^. This information will provide a high- 
precision determination of the expansion history of the Uni- 
verse through accurate measurements of the BAO feature in 
the large-scale galaxy clustering. BOSS will also attempt to 
obtain, for the first time, BAO measurements at high red- 
shift {z « 2.5) through the Lya forest absorption spectra of 
about 150,000 quasars. 

The increasing precision of the new surveys demands 
accurate models of the LSS observations to extract the maxi- 
mum amount of information from the data without introduc- 
ing biases or systematic effects. The BAO signal in the corre- 
lation function and power spectrum is modified by the non- 
linear evolution of dens ity fluctuations, redshift-space distor- 
tions , and galaxy b ias ('Mci ksin et al.lll999l:lEisenstein et al 



CMASS sample f Anderson ot al."2012': 'Ma nera et~ai] 12012 



m prer 



20071 : ISeo fc EiseiTs tcin 20071: ISeo et al.l I2OO8I: 



lAngulo et all 12003: ICrocce 



Sanchez et al.l I2OO8I: [Cott et all [2009 : Kim et al] l2009l : 
Montesano et al.ll20ld : lKim et al1l201ll ). These effects must 



iSmith et aL. 

Scoccimarroj 20081: 



be taken into account in the models used to interpret 
the observations. New developments in perturbation the- 
ory, such as Renormalized Perturbation Theory (RPT, 
ICrocce fc Scoccimarrd 120061 ) . have provided substantial 
progress regarding the theoretical understanding of the 
effects of non-linear evolution, which can now be accu- 
rately modelled IC roccc &: Scoccimarro 2006; Matsubara 
l2008al bl: Taruva et al. 2009) , and even partially corrected for 
(lEisenstein et al.ll2007. : .Seo et al. 20ld: Padmanabhan et al.l 
|2012| ). Based on RPT. ICrocce fc Scoccimarrd (|2008l ) pro- 
posed a simple model to describe the full shape of the 
correlation function on large scales. [Sanchez etal.| (|2008l ) 
showed that this model yields an excellent description 
of the results of N-body simulations, providing a ro- 
bust tool to extract unbia sed cosmologic a l cons traints out 
of measurements of £,(r). I Sanchez et al.l l|2009l ) used this 
model to obtain constraints on cosmological parameters 
from the co rrelation function of a sample of LRGs from 
SD SS-DR6 (lAdelman-McCarthv et al.l I2OO8I ) as measured 
bv ICabre fc Gaztanaga 1120091 ). The same ansatz h as been 
used bv lBeutler et all l|201ll ) and lBlake et af] l|20lj ) for the 
analysis of the correlation functions of the 6dF and Wig- 
gleZ galajcy surveys, respectively. A n analogous approach 
was used bv lMontesano et al.l (|2012l ) to study the cosmolog- 
ical implications of the LRG power spectrum in SDSS-DR7 
l|Abazaiian et al.|[2009l ). 

In this paper we ap ply the parametrization of 
ICrocce fc Scoccimarrd (|2008l ) to the redshift-space correla- 
tion function of a high redshift galaxy sample from BOSS 
Data Release 9 (DR9). This sample, denoted CMASS, is con- 
structed through a set of colour-magnitude cuts designed 
to select a rough ly volume-limited sam ple of massive, lu- 
minous galaxies IjEisenstein et al.l I2OIII . Padmanabhan et 
al. in prep.). We combine the CMASS clustering informa- 
tion with recent measurements of CMB, BAO and type 
la supernovae data. We derive constraints on the param- 
eters of the standard ACDM model, and on a number of 
potential extensions, with an emphasis on the constraints 
on the dark energy equation of state. Our analysis is part 
of a series of papers aimed at providing a thorough and 
comprehensive description of the galaxy clustering in the 



Rcid ct al. 2012: Ross ct al. 2012 : Samushia et all 
iToieiro et al.ll20li )! 

The outline of this paper is as follows. In Section [5] we 
describe the galaxy sample that we use and the procedure we 
follow to compute its correlation function. We also present a 
discussion on the cosmological information contained in this 
measurement. Section [3] describes the additional datasets 
that we combine with the CMASS correlation function to 
obtain constraints on cosmological parameters. Our model 
of the full shape of the correlation function, the parame- 
ter spaces we explore and the applied methodology is de- 
scribed in Section [31 In Section [S] we present our results for 
constraints on cosmological parameters from different com- 
binations of datasets and parameter spaces. In Section [6] 
we analyse the differences in the clustering of the northern 
and southern Galactic hemispheres and explore their im- 
plications on the obtained cosmological constraints. Finally, 
Section [7] contains our main conclusions. 



2 CLUSTERING ANALYSIS OF THE 
BOSS-CMASS GALAXIES 

We base our analysis on the large-scale two-point correla- 
tion function, ^(s), of the BOSS-CMASS galaxy sample. In 
this Section we review the most important details of the 
construction of the sample (Section [2TT}, and our clustering 
analysis fSection l2.2p . 



2.1 The CMASS galaxy sample 

The galaxy target selection of BOSS consists of two sepa- 
rate samples, dubbed LO WZ and CMASS, desig ned to cover 
different redshift ranges (jEisenstein et al.ll201ll . Padmanab- 
han et al. in prep.). These samples are selected on the basis 
of photometric obse rvations done wit h the dedicated 2.5- 
m Sloan Telescope (jGunn et al.l 120061 ). located at Apache 
Point Observatory in New Mexico, us ing a drift-scanning 
mosaic CCD camera (jGunn et al.lll998l ). These samples are 
constructed on the basis of gri colour cuts designed to select 
luminous galaxies at different redshifts at a roughly constant 
number density. Spectra of the LOWZ and CMASS samples 
are obtained using the double-armed BOSS spectrographs, 
which are significantly upgraded from those used by SDSS- 
I/II, covering the wavelength ra nge 3600 A to 10000 A with a 
resolving power of 1500 to 2600 (|Smee et akllin prep.l ). Spec- 
troscopic redshifts are then measured u sing the minimum- 



i mum-Y 

template- fitting procedure described in lAihara et al.l (|201ll ). 
with template s and methods upd ated for BOSS data as de- 
scribed in Bol ton et al.l din prep.!) . 

Our analysis is based on the clustering properties of the 
CMASS sample, which is selected to be an approximately 
complete galaxy sample down to a limiting stellar mass 
fMa raston et al.l lin prep.!) . The CMASS sample is domi- 
nated by early type galaxies, although it con tains a signifi- 
cant fraction of massive spirals (~26 per cent. [Masters et al.l 
I2OIH ) . Most of the galaxies in this sampl e are central galax- 
ies, with a ~10 p el 
iNuza et al.ll2012l). 



20 111 ). Most ot the galaxies m this sampl e are central galax- 
ies, with a ~10 p er cent satellite fraction l| White et al.ll201ll : 



[Anderson et al.l l|2012l ) presents a detailed description 
of the construction of the catalogue for LSS studies based 
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Completeness 
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0.85 0.90 0.95 1.00 





Figure 1. The sky coverage, in Galactic coordinates, of the CMASS DR9 spectroscopic sample used in this analysis in the northern 
(left panel) and southern (right panel) Galactic hemispheres. Different sectors are colour-coded according to their completeness. The 
low completeness at many edges is due to planned but currently unobserved tiles that will overlap the current geometry. The light grey 
shaded region shows the expected footprint of the final survey, totalling 10,269 deg^ . 



on this sample, and the calculation of the completeness of 
each sector of the survey mask, that is, the a reas of the sky 
cover ed by a unique set of spectroscopic tiles IjBlanton et alj 
12003). which we characterize usinR the Mangle software 
l|Hamilton fc Tegmark|[20o3 : ISwanson et al.ll2008l ). We only 
include sectors with completeness larger than 75%. Our re- 
sults are not affected by this limit, as this leaves out only a 
small fraction of the total DR9 area. We restrict our anal- 
ysis to the redshift range 0.43 < z < 0.7, producing a final 
sample of 262,104 galaxies, of which 205,947 and 56,157 are 
located in the Northern and Southern Galactic hemispheres, 
respectively. Fig. [T] shows the angular footprint, in Galac- 
tic coordinates, of the resulting sample for the Northern 
(left) and Southern (right) Galactic caps (hereafter NGC 
and SGC, respectively), colour-coded according to sector 

com pleteness. 

iNuza et al.l (|2012l ) compared the small and 
intermediate-scale clustering of this sample to the ex- 
pectations of a flat ACDM cosmological model by applying 
an abundance matching technique to the Multida r k sim - 
ulation. In thr e e comp anion papers, Reid et al.l l|2012h . 



ISamushia et al.l (|in predl and iToieiro et al.l (|2012l 'l study 



the signature of redshift-space distortions in this sample 
and explore its cosmological implications. Here we focus on 
the shape of the large-scale monopole correlation function 
to obtain constraints on cosmological parameters. 



2.2 The redshift-space correlation function 

We characterize the clustering of the CMASS galaxy sam- 
ple by means of the angle-averaged redshift-space two-point 



correlation function £,{s). Here we summarize the procedure 
we follow to obtain this measurement. 

The first step in the calculation of three-dimensional 
clustering statistics is the conversion of the observed red- 
shifts into distances. For this we assume a fiat ACDM fidu- 
cial cosmology with matter density, in units of the critical 
density, of 0.-^ ~ 0.274, and a Hubble parameter h — 0.7 (ex- 
pressed in units of 100 km s~^Mpc~^). This is the same fidu- 
cial cosmology a ssumed by White et al.l (12011 ) and our com- 
panion p apers 'Anderson et al. 201Z; Manera et al.| |2012| : 
I Ross et al . 2012; Reid ct al. 2012; Tojciro et al. 20li). As 
will be discussed in Section 14.31 the choice of the fiducial 
cosmology has implications on the resulting correlation func- 
tion. 

We then compute the full correlation function ^{s,fi), 
where fi = S||/|s| and S|| is the radial c omponent of the 
separa tion vector s, using the estimator of iLandv fc Szalavl 
(Il993h . namely 



C(s,m) = 



DP - 2DR + RR 
RR ' 



(1) 



where DD, DR and RR are the normalized pair counts in 
each bin of (s, n) in the data and a random sample with 50 
times more objects than the original data, constructed to 
follow the same selection function (for more details on the 
const ruction of the random catalogue see Anderso n et al.l 
[2013). We infer the angle- averaged redshift-space correlation 
function as the monopole of (,{s,fi), that is 



(2) 



C(s) = 2 / ?(«'^)dM- 



This method should be preferred over the commonly used 
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one, in which the DD, DR, and RR counts are integrated 
over /i before they are combined as in equation ([l| to com- 
pute £,(s), ignoring the fact that the geometry of the surve; 
introduces a /i d ependence on RR ( Samushia et alj 1201 



I 



iKazin et al.ll2012l ) , although the differences between the two 
approaches are more significant for higher multipoles. 

When computing the pair counts in equation ([T]), a few 
important corrections must be taken into account. This is 
done by assigning a series of weights to each object in the 
real and random catalogues. First, we apply a radial weight 
given by 

"- ^ 1 + Pln(z) ' (') 

where n{z) is the expected number density of the cata- 
logue at t he gi ven redshift and P-u, is a free parameter. 
iHamiltonl (| 19931 ) showed that setting P^ = A-kJ-i,{s), where 
Ji{s) = /q ^(s')s''^ds', minimizes the variance on the mea- 
sured correlation function for the given scale s. Follow- 
ing standard practice we use a scale-ind ependent value of 
P™ = 2 X 10"* fe~^Mpc^. lReid et al.1 (|2012l ) show that the full 
scale-dependent weight provides only a marginal improve- 
ment over the results obtained using this constant value. 

We include additional weights to account for non- 
random contributions to the sample incompleteness and to 
correct for systematic effects. The incompleteness in a given 
sector of the mask has a random component due to the fact 
that not all galaxies satisfying the CMASS selection criteria 
are observed spectroscopically. In any clustering measure- 
ment this is taken into account by down-sampling the ran- 
dom catalogue in that region of the sky by the same fraction. 
However, there are two other sources of missing redshifts 
which require special treatment: redshift failures and fibre 
collisions. 

Even when the spectrum of a galaxy is observed, it 
might not be possible to obtain a reliable estimation of the 
redshift of the obje ct, leading t o wh at is called a redshift 
failure. As shown in I Ross et al.l l|2012l ). the probability that 
a spectroscopic observation leads to a redshift failure is not 
uniform across the field since these tend to happen for fibres 
located near the edges of the observed plates. Hence, these 
missing redshifts cannot be considered as an extra compo- 
nent affecting the overall completeness of the sector. 

However, the main cause of missing redshift is fibre col- 
lisions (Zehavi et al. 2002 ; Masiedi et al. 2006 ). The BOSS 
spectrographs are fed by optical fibres plugged on plates, 
which must be separated by at least 62". It is then not 
possible to obtain spectra of all galaxies with neighbours 
closer than this angular distance in one single observation. 
The problem is alleviated in sectors covered by multiple ex- 
posures but, in general, it is impossible to observe all the 
objects in crowded regions. 

To correct for these effects we follow IRoss et al.l (|2012l ) 
and implement two sets of weights, w^i and wtc, whose de- 
fault value is 1 for all galaxies in the sample. For every galaxy 
with a redshift failure, we increase by one the value of w-d of 
the nearest galaxy with a good redshift measurement. Sim- 
ilarly, for each galaxy whose redshift was not observed due 
to fibre collisions, the value of Wic of its neighbour, closer 
than 62" , is increased by one. These are then combined into 
a single weight to correct for missing redshifts given by 
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Figure 2. Panel (a): spherically averaged redshift-space two- 
point correlation function of the full CMASS sample. The error- 
bars were obtained from a set of 600 mock catalogu es constructed 
to fol low the same selection function of the survey JManera et al.l 
l2012t) . The dashed line corresponds to the best- fitting ACDM 
model obtained by combining the information from the shape 
of the correlation function and CMB measurements (see Section 
15.11 1. Panel (b): same as panel (a), but rescaled by (s/sbao)^i 
where sbao = 153.2 Mpc (which corresponds to 107.2 h~^}Apc), 
to highlight the baryonic acoustic feature. 



On the scales analysed in this paper, the application of these 
weights effectively corrects for the effects of fibre collisions 
and redshift failures providing an excellent agreement with 
the results obtained using the method recently proposed by 
IGuo et al.l (l201lh. 

IRoss et al.l ([2013) performed a detailed analysis of the 
systematic effects that could potentially affect any clustering 
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measurement based on the CMASS sample showing that, 
besides redshift failures and fibre colhsions, other important 
corrections must be considered in order to obtain unbiased 
clustering measurements. They found that the local stellar 
density is the dominant source of systematic errors as it has 
a significant effect on the probability of detecting a CMASS 
galaxy. In this way, the variations of stellar density across 
the sky introduce spurious fiuctuations in the gal axy density 
field w hich affect all clustering measurements. IRoss et al.l 
lJ2012t ) found that this systematic effect can be corrected 
for by applying a set of weights w^ys which depend on both 
the stellar density and the galaxy iflbcr2 magnitude, that is, 
the i-band magnitude measured within a 2" aperture. We 
include these weights in the final total weight lutot used in 
all our clustering measurements 



^'' " {N^ - 1) 



E(a(si)-C>i))(a(sj)"e>j)), (6) 



Wtot = lUrlUmr Ways. 



(5) 



Additional potential systematics such as Galactic ex- 
tinction, seeing, airmass, and sky background have also 
been investigated, and all have been found to potentially 
introduce much smaller spurious fluctuations. These non- 
cosmological fluctuations can be corrected for using a 
weighting scheme that minimises these fluctuations as a 
function of a given systematic effect. 

The upper panel of Fig. [2]shows the large-scale redshift- 
space correlation function of the full CMASS sample ob- 
tained through the procedure described above. The dashed 
line corresponds to the best-fitting ACDM model obtained 
from the combination of this measurement with CMB ob- 
servations as described in Section 15.11 The BAO peak can 
be seen more clearly in the lower panel, which shows the 
same measurement rescaled by the ratio (s/sbao)^, where 
sbao = 153.2 Mpc corresponds to the sound horizon scale 
in our fiducial cosmology. As will be discussed in more de- 
tail in Section [6l the measurements of the two-point corre- 
lation function in the NGC and SGC sub-samples exhibit 
intriguing differences. Although the overall shapes of these 
measurements are similar, they show differences at the scale 
of the acoustic peak. In Section |6] we discuss the significance 
of these differences and their impact on the inferred cosmo- 
logical constraints. 

To obtain an estimate of the covariance matrix of the 
correlation function measured in these samples, we use a 
set of A'm = 600 independent mock catalogues based on a 
method similar to PTHalos ('Scoccimarro fc Shethl l2002l l. 
although with some important differences. A detailed de- 
scription of the construction of these mock catalogues and 
a compari son with the r e sults o f N-body simulations is pre- 
sented in iManera et al.l (|2012l l^l. These simulations corre- 
spond to the same fiducial cosmology used to measure ^(s) 
in the real catalogue and were designed to follow the selec- 
tion function of the NGC and SGC CMASS sub-samples. 
We measured the correlation function of each mock cata- 
logue using the same binning schemes as for the real data 
and the radial weights of equation ([3]). We then use these 
measurements to obtain an estimate of the covariance ma- 
trix of i{s) in the NGC, SGC as 



^ These mock catalogues will be made available in 
|http://www.marcmanera. net/mocks/I 



where ^k(s) is the correlation function from the fe-th mock 
catalogue, and ^(s) is the mean correlation corre lation func- 
tion f rom the ensemble of realizations. As in IRoss et al.l 
l|2012f ). we assume that the NGC and SGC regions are in- 
dependent and compute the covariance matrix of the full 
CMASS sample as Cj^J^ — C^qq + ^SGC- T'^® errorbars 
in Fig. [2] correspond to the square root of the diagonal en- 
tries in Cfuii. 



3 ADDITIONAL DATA-SETS 



As described in Section [4.31 the two-point correlation func- 
tion contains valuable cosmological information. However, 
it is not possible to constrain high-dimensional parameter 
spaces to high precision using this measurement alone. Here 
we describe the additional datasets with which we combine 
the CMASS ^(s) in order to improve the obtained cosmo- 
logical constraints. 

Undoubtedly, the measurements of the temperature and 
polarization fluctuations of the CMB constitute the most 
powerful and robust cosmological probe to date. In partic- 
ular, the results f rom the seven-year o f observations of the 
WMAP satellite ('Hinshaw et al."2009') and the South Pole 
Telescope (SPT, Keisler et al. 2011) provide a detailed pic- 
ture of the structure of the acoustic peaks in the CMB power 
spectrum up to multipoles I ~ 3000. This information places 
tight constraints on the parameters of the basic ACDM 
model. However, the power of these observations is limited 
by nearly exact degeneracies that arise when deviations from 
this simple model are explored (JEfstathiou fc Bondlll999l ). 
These degeneracies can be broken by combining the CMB 
information with additional datasets, such as the shape of 

In our analysis we use the temperature power spec- 
trum in the range 2 ^ f ^ 1000 and the temperature- 
polarization power spectrum for 2 ^ ^ ^ 450 fr om the seven- 
year of observations of the WMAP satellite ( Jarosik et al.l 
I2OIII : iKomatsu et al.l I2OIII : iLarson et al] boill). combined 
with the recent SPT observations of iKeisler et al.l (12011! ) 
for 650 < f < 3000. While for i < 650 the CMB power 
spectrum is dominated by primary anisotropics, at smaller 
angular scales it contains a non-negligible contribution from 
secondary anisot ropies. To take this i nto account, we follow 
the treatment of iKeisler et al.l l|201H ) and include the con- 
tributions from the Sunyaev-Zel'dovich (SZ) effect, and the 
emission from foreground galaxies (considering both a clus- 
tered and a Poisson point source contribution) in the form of 
templates whose amplitudes are considered as nuisance pa- 
rameters and marginalized over. These templates are only 
applied to the SPT data. We refer to the WMAP-SPT com- 
bination as our "CMB" dataset. 

Additionally, we consider the constraints provided by 
the Hubble diagram of type la super n ovae (SN) obtained 
from the compilation of IConlev et al.l (|201ir ). This sample 
contains 472 SN, combining the high-redshift SN from the 
flrst three years of the Supernova Legacy Survey (SNLS) 
with other samples, primarily at lower redshifts. In order to 
take into account the effect of the systematic errors in our 
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cosniological constraints we follow the recipe of lConlev et alj 
l|201]| |. who performed a detailed analysis of all identified 
systematic uncertainties, characterizing them in terms of a 
covariance matrix that incorporates effects such as the re- 
cently discovered correlations between SN luminosity and 
host galaxy properties, as well as the uncertainties of the 
empirical light-curve models. When only SN data are used to 
constrain cosmological parameters, the uncertainty budget 
is dominated by statistical errors. However, when these data 
are combined with external datasets, as in our case, statisti- 
cal and systematic uncertainties are comparable, highlight- 
ing the importance of an accurate treatment of the later. 

We also use information from other clustering mea- 
surements in the form of constraints on i/s(z) and A{z) 
fro iu independen t BAO analyses. We use the results 
of iBeutler et al.l (|201lh . who obtained an estimate of 
ys{z = 0.106) = 0.336 ± 0.015 from the large-scale 
correlation funct ion of the 6dF Galaxy Survey (6DFGS, 
I Jones et al.ll2009l ). We also include the 2% distance mea- 
surement of (i/B(0.35) ) ~^ — 8.88 i: 0.17 recently o btained by 
iPadmanabhari et al.l (|2012l ) and IXu et all (|2012l ) from the 
application of an u pdated version of the re construction tech- 
nique proposed bv lEisenstein et ahl l|2007l 'l to the clustering 
of galaxies from the final SDSS-II LRG sample. The applica- 
tion of this algorithm resulted in an improvement of almost 
a factor two in the accuracy on ys over the constraint ob- 
tained from the unreconstructed sample. We combine the 
result from th ese analyses into o ur "BAO" dataset. In a re- 
cent analysis, iBlake et alj |201lh used the full shape of the 
two-point correlation function f rom the final dataset of the 
WiggleZ Dark Energy Survey IjDrinkwater et al.l 120101 ) to 
obtain constraints on ys{z) and A{z) for three independent 
redshift slices of width Az = 0.4. We do not include these 
measurements in our analysis given the significant overlap 
of the WiggleZ data with the sample a nalysed here. How- 
ever, as shown in [Anderson et al.l (|2012l ). the WiggleZ BAO 
measurements are in excellent agreement with those inferred 
from the CMASS sample. 

The datasets described above are used in different com- 
binations to check the consistency of the constraints re- 
turned. We start from the constraints obtained using CMB 
data alone, which we then combine with the CMASS cor- 
relation function in our "CMB-I-CMASS" combination. We 
then add separately the SN and additional BAO data to 
test the impact of these datasets on the obtained results. 
Our tightest constraints are obtained from the combination 
of all four datasets. 



4 METHODOLOGY 

We obtain constraints on cosmol ogical paramet e rs foU ow- 
ing a similar approach to that of I Sanchez et al.l ((20091). In 
this Section we summarize the main points of our analysis 
method. The parametric model we use to describe the shape 
of the correlation function in redshift space is summarized 
in Section 14.11 Section 14.21 describes the different parame- 
ter sets that we consider, together with the methodology we 
follow to explore and constrain them. Section [4.31 describes 
the way in which cosmological information is extracted out 
of a measurement of ^(s). 



4.1 Modelling the full-shape of ^(s) 

Following i Crocce fc Scoccim arrd (|2008l ) and I Sanchez et al.l 
(2008), we model the shape of the large-scale correlation 
function, (,{s), by applying the following parametrization: 

C(s) = 6' [Us)S)e-^'*^^"+AMC^Us)^i'Hs)\, (7) 

where b, fc* and Amc are treated as free parameters, and the 
symbol (g) denotes a convolution. Here ^l is the derivative of 



the linear correlation function ^l, and ^^ is defined by 



(1) 



^«(s)EES.V-leL(s) 



27r2 



Phik)ji{ks)kdk, 



(8) 



with ji{y) denoting the spherical Bessel function of or- 
der one. This parametrizatio n was originally proposed by 
ICrocce fc Scoccimarrd (|2008l) and it is based on the the- 
oretical &;am«work_of_^£iwrmalized Perturbation Theory 
(RPT, ICrocce fc Scoccimarrd 120061 ) . where the non-linear 
power spectrum Pj-ii,{k,z) can be computed as the sum of 
two terms 



PNL(fc, z) = G{k, zfPl^ik, Z) -f PMc(fc, Z). 



(9) 



The first of these contributions represents a re-summation 
in the renormalized propagator, G(fe, 2), of all the terms in 
the perturbation theory expansion of PNL(fc, z) proportional 
to the linear theory power spectrum P\,{k\ The second term 
groups all the remaining contributions, which arise from the 
coupling of different Fourier modes. The non-linear correla- 
tion function is then given by an analogous decomposition, 
which motivates the parametrization of equation ((Tj. The 
exponential in the first term of equation ([7]) is based on the 
fact that, in the high-fc limit, the propagator can be accu- 
rately described as a Gaussian damping, while the second 
term corresponds to the leading order contribution to ^mc 
arisi ng from the co upling of two initial modes. 

ISanchez et al.l (^2008) compared this model against the 
results of an er isemble of large volume N-body simulations 
('L-BASICC-H. lAngulo et al.ll2008i ') at various redshifts, and 
showed that it provides an accurate description of the full 
shape of the correlation function, inclu ding also the effects 
of bias and redshift-space distortions. Sa nchez et al.l (|20091 ) 
used this model to obtain constraints on cosmological pa- 
rameters from the correlation function of the LRG sampl e 
from SDSS-DR6 measured by ICabre fc Gaztanagal (|2009|) 
This parametrization has also been used by iBeutler et al.l 
(2011) and Blake ct al. (2011) for their analyses of the cor- 
relation function measurements from the 6dF Galaxy Sur- 
vey an d the WiggleZ Dark Energy Survey. iMontesano et al.l 
(|2012t) applied an analogous parametrization to study the 
cosmological implications of the power spectrum of an LRG 
sample drawn from SDSS-DR7. 

The smoothing length fe* depends on cosmology and 
redshift, but also on galaxy type through its dependence on 
halo mass. For this reason, we follow a conservative approach 
and consider kj, as a free par ameter. 

FoUowing lSanchez et al.l ((2009), we restrict the compar- 
ison of the model of equation ((T} and the measured CMASS 
correlation function to s > 40/i~^Mpc. Although this is a 
conservative lower limit, on smaller scales further contri- 
butions to ^Mc(s) should be considered. We also limit our 
analysis to scales s < 200/i~^Mpc, since on larger scales 
all viable models predict similar shapes for ^(s). We com- 
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pute the likelihood of the model assuming a Gaussian form 
C PC exp( — y^/2). T his choice is justified by the results of 
iManera et alj l|2012h . who found that the probability distri- 
bution function of ^(s) inferred from the ensemble of mock 
catalogues can be described by a Gaussian distribution to 
high accuracy. 

To allow for the fact that, when computing the CMASS 
correlation function, galaxy distances were calculated with 
our fiducial cosmology, a correction must be applied to the 
model before computing its corresponding y^ value (see Sec- 
tion [331). 



4.2 Cosmological parameter spaces 

The starting point of our analysis is the basic ACDM param- 
eter space. This is the simplest model able to successfully 
describe a large variety of cosmological datasets. It corre- 
sponds to a flat universe where the energy budget contains 
contributions from cold dark matter (CDM), baryons, and 
dark energy, which is given by vacuum energy or a cosmo- 
logical constant A (i.e. with an equation of state parameter 
u;de = — 1). Primordial density fluctuations are adiabatic, 
Gaussian, and have a power-law spectrum of Fourier am- 
plitudes, with a negligible contribution from tensor modes. 
This model can then be defined by specifying the values of 
the following six parameters: 



PaCDM = (l^b,tt^dm, 6,T", As 



(10) 



The baryon and dark matter densities, Ub = 0.hh^ and 
Wdm = fidm^^, and the ratio between the horizon scale at 
recombination and the angular diameter distance to the cor- 
responding redshift, 0, characterize the homogeneous back- 
ground model. This set is equivalent to fixing the values of 
f^b, fidm and h, but it is better constrained by the CMB 
data. The primordial power spectrum of the scalar fluctua- 
tions is described by its amplitude. As, and spectral index, 
Wb. The values of these parameters are quoted at the pivot 
wavenumber of A; = 0.05 Mpc~^. Finally, r gives the optical 
depth to the last scattering surface, which we compute as- 
suming instantaneous reionisation. Our constraints on the 
ACDM parameter space are described in Section FS. II 

In order to constrain possible deviations from the 
ACDM model, in Sections 15. 2145. 51 we explore a number of 
extensions of this parameter space by allowing for variations 
on the following set of parameters: 



Pcxtra = (f2k,/i/,r, Wde). 



(11) 



These are the curvature of the Universe, the dark matter 
fraction in the form of massive neutrinos, /^ = Sliz/Sldm, 
the tensor-to-scalar mode amplitude ratio of the primordial 
fluctuation^, and the dark energy equation of state param- 
eter. For most of this paper, we assume that the dark energy 
equation of state is independent of redshift. In Section [5.5.21 
we allow also for a time variation of this parameter using 
the st andar d linear param etrization of lChevallier fc Polarski 
l|200ll ) and lLindeiJ J2003l ) given by 



lODE(a) = Wo + Wa(l — a). 



(12) 



^ When including tensor modes we assume the slow-roll consis- 
tency relation and fix the tensor spectral index as nt = — r/8. 



where a is the expansion factor and wo and Wa are the pa- 
rameters we constrain. 

We also present constraints on other quantities which 
can be derived from the ones listed in equations (|10p and 
ITT}. In particular we are interested in: 



(SIde, f^m, h,as,to, Zrc , -Pv (^m) , /^ m„ , /) . 



(13) 



These are the dark energy and total matter densities (i.e., 
including the contributions from baryons, cold dark matter 
and neutrinos), the Hubble factor, the rms linear pertur- 
bation theory variance in spheres of radius 8/i~^Mpc, the 
age of the universe, the redshift of reionization, the aver- 
age distance to the mean redshift of the sample (given by 
equation I17|l , the sum of the neutrino masses, given by 



^ m^ = 94.4 cjdm/^ eV, 



(14) 



and the logarithmic derivative of the growth factor, f{zm) ~ 
dlnD/dlna. 

We explore these paramete r spac es using the Cos- 
MOMC code of iLewis fc Br"idi3 (|2002| ). COSMOMC uses 
CAMB to co mpute power spec tra for the CMB and matter 
fluctuations (JLewis et al.l200d ) . We use a generalized version 
of CAMB which supports a time-dependent dark energy equa- 
tion of state (.F ang et al.n200&). We incl uded additional mod - 
ifications from lKeisler et al.l ()201ll l and lConlev et al.l (|2011h 
to compute the likelihood of the SPT and SNLS datasets. 

In order to compare a given cosmological model with 
the datasets described in Sections [2] and [3] it is necessary to 
include a set of nuisance parameters given by 



nuisance 



(6, fc*, AmC,-D3000,J 



, -C>3ooo,a,/3), 



(15) 



to the parameter sets described above. The bias factor b, the 
damping scale fc* and the mode-coupling amplitude Amc 



are described in detail in Section |4TT] The quantities -Dtoooi 
^3000 and -D3000 gi"^6 Ihs amplitudes of the contributions 
from the Sunyaev-Zel'dovich effect, the clustering of the 
foreground emissive galaxies and their shot-noise fluctua- 
tion power, respectively, to the high-^ CMB angular power 
spectrum. The foreground terms are used only when calcu- 
lating the SPT likelihood; they are not used when ca l culat - 
ing the WMAP likelihood. We follow iKeisler et all (I2OIII ) 
and apply Gaussian priors on the amplitude of each of 
these foreground terms given by -D3(foo ~ 19.3 ± 3.5 fiK^, 
DsMO = 5.0 ± 2.5 pK^ and L»f|oo = 5.5 ± 3.0 /xK^ The 
parameters a and /3 are additio nal nuisance parameters in- 
troduced bv lConlev et al.l (120111 ) for the correct treatment of 
the systematics in the analysis of the SN data. When quot- 
ing constraints on the parameters of equations (|10p -(|li p . the 
values of these parameters are marginalized over. 



4.3 Extracting information out of ^(s) 

In this Section we describe the information encoded in the 
shape of the two-point correlation function and how it can 
be used to obtain constraints on cosmological parameters. 
As described in Section [2]2l the measurement of the correla- 
tion function requires the assumption of a fiducial cosmology 
to map the observed redshifts into distances. This fact has 
important implications on the parameter combinations that 
are constrained by £,{s). 

Different choices of the fiducial cosmology lead to a 
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Figure 3. Panel (a): mean correlation function from our en- 
semble of mock catalogues obtained by assuming the true cos- 
mological parameters as fiducial values (circles connected by a 
solid line) and a flat cosmology with Hm = 0.4 (squares con- 
nected by a dashed line). The shaded region correspond to the 
variance between the dilferent realizations of the ensemble. Panel 
(b); same measurements as panel (a), but expressed as a function 
of 2/ = s/D-v{zm), which removes the dependence on the fiducial 
cosmology. 



rescaling of the distances to the individual galaxies s — >■ 
s' , affecting the volume element of the survey. This effect 
can be encapsulated in the Jacobian of the transforma- 
tion llEisenstein et al.ll2005l : ISanchez et al.ll2009l : iKazin et al.l 
|2012| ) 



d^s' = (^^4H ] d's 
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(16) 



Here Dv{zni) is a measure of the average distance to the 
mean redshift of the survey, «m = 0.57, given by 



Dviz) = {{i + zyDA(zy 



H{z) 



1/3 



(17) 



where H{z) is the Hubble parameter and Da (2) is the proper 
angular diameter distance. 

Fig. |3] illustrates the effect of assuming different fidu- 
cial cosmologies on the measurement of f(s). The points 
connected by a solid line in panel (a) show the mean corre- 
lation function of our ensemble of mock catalogues, obtained 
assuming as fiducial cosmology the true values of the simula- 
tion parameters. The shaded region corresponds to the vari- 
ance between the individual realizations. The squares con- 
nected by a dashed line correspond to the mean correlation 
function from the same set of mock catalogues, but obtained 
assuming a fiat ACDM model with Q.^ = 0.4. The two mea- 
surements show significantly different slopes and positions 
of the acoustic peak. As equation (|16|) suggests, this change 
is simply due to a rescaling of the horizontal axis. This effect 
can be better appreciated in panel (b) of Fig. [S] where the 
impact of the fiducial cosmology has been removed by ex- 
pressing the measured correlation functions in terms of the 
dimensionless variable y = s/Dv'^(zni). This exercise shows 
that, although the true underlying correlation function is 
not a real observable, it is possible to obtain a measurement 
which is independent of the fiducial cosmology by expressing 
it as i{y). 

The particular choice of the fiducial cosmology must be 
taken into account when comparing a measurement of ^(s) 
with theoretical predictions. As described above, this can 
be achieved by expressing both model and measurements in 
terms of y. Alternatively, the effect of the fiducial cosmology 
might be introduced in the model by rescaling the scales s 
by a factor 



D^''""'"(2m) 
r)niodcl 



(Zm) 



(18) 



before comparing it to the measured ^(s). We follow this 
approach in our analysis. 

The most important source of cosmological information 
in ^(s) is the location of the acoustic peak, which is closely 
related to the sound horizon at the drag redshift rs(2d). As- 
sociating the position of the peak in ^(j/) with this scale, it 
is clear that this measurement will provide constraints on 
the parameter combination 



ys{Zr, 



^{Zd) 



Dv{z 



(19) 



However, the location of the acoustic peak does not cor- 
respond exactly to the acoustic scale. Non-linear evolu- 
tion and redshift space distortions damp the a coustic peak 
and sh ift its position toward s smal le r scales |Smith et al.l 



2008 



1 ICrocce fc Scoccimarrol I2OO8I : lAngulo et all I2OO8I : 



'San chez et al.ll2008l ). Nonetheless, if these effects are mod- 
elled correctly, a measurement of ^(s) would still provide 
constraints on the parameter combination of equation (|19p . 
allowing for the correct underlying cosmology to be recov- 
ered. 

Similarly, the measurement of the power spectrum, 
P{k), will be subject to the same effect, which can be 
removed by multiplying the measured wave-numbers by 
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Figure 4. The 68 and 95 per cent marginalized constraints in 
the i^ni — Dy{zui) plane, where Wm = Slm^i^, obtained from the 
shape of the CMASS correlation function alone (solid lines). The 
dashed, dot-dashed and dotted lines correspond to constant values 
oi Dv{zra)^m, ys{2m) (equation I19I I. and A{zin) (equation I20I I. 
respectively. 



£'v(zm)- In this way, the wavelength of the acoustic oscilla- 
tions inferred from a measurement of P{k) will provide con- 
straints on 2/7^(2m)- As P{k) is not a dimensionless quan- 
tity, its amplitude is also affected by the fiducial cosmol- 
ogy (by a factor proportional to (73v(zm))^). This can be 
avoided by working with the dimensionless power spectrum 
A2(fc) = P(fc)fcV(27r2). 

Besides the BAO, the power spectrum contains informa- 
tion on another useful scale. The location of the turn-over in 
P(k) is related to the size of the sound horizon at the time of 
matter-radiation equality. In the absence of massive neutri- 
nos, and for a fixed effective number of relativistic species, 
this scale is fceq oc fim^'^Mpc"^. Taking into account the 
effect of the fiducial cosmology, the quantity that can ac- 
tually be constrained is fccqDv(zm). The information about 
this parameter combination is also encoded in the shape of 
the correlation function, where it is related to the position 
of the zer o-crossing at scale s larger than those of the acous- 
tic peak (jPrada et al.ll201lh . In this way, a measurement of 
^(s) provides constraints on the same parameter combina- 
tion. This quantity is degenerate with other parameters, like 
the baryon density and the scalar spectral index, which also 
affect the shape of ^(s). However, the later are tightly con- 
strained by CMB observations (e.g.. iKomatsu et al.ll201ll : 
iKeisler et a"i]|201ll ). 

The contours in Fig. [4] show the two-dimensional 
marginalized constraints in the cJm — Dv(zm) plane, where 
cdm = 0.mh? ^ obtained from the shape of the CMASS corre- 
lation function, using the model described in Section[4T] To 
ameliorate the effect of the degeneracies between uim, and cjb 
and Us in this exercise, we have applied Gaussian priors of 
Ub = 0.0222±0.0010 and n^ = 0.966±0.020. These priors are 



weaker than the corresponding accuracy with which these 
parameters are determined by current CMB data (see Sec- 
tion [5]), allowing us to quantify more clearly the information 
provided by ^(s). The full combination of this measurement 
with CMB data will result in slightly tighter constraints. 

The dashed and dotted lines in Fig. |4] correspond to 
constant values of i/s(zin) and Dv(zn-i)ctJn-i- The interplay 
between the constraints on these parameter combinations 
shapes the allowed region in the cum — Dy plane. The dotted 
line in Fig. [4] correspon ds to a constant value of the quantity 
dEisenstein et al.ll2005l 'l 



A{Zm) = Dv(Zm) 



yf^mJ^o 



(20) 



which approximately describes the resulting degeneracy be- 
tween LJm and Dv(zm). To take into account the effect of the 
baryon density on the scale fccq, this quantity should be de- 
fined in terms of the shape parameter P. However, we main- 
tain the usual definition to simplify the comparison with 
previous analyses. The CMASS correlation function implies 

a constraint of A{zm) = 0.444 ± 0.014 

Two of ou r comp anion papers, Ander son et al.l (|2012l ) 
and lReid et al.l l|2012r ). study the cosmological implications 
of the galaxy cluster ing in the CMASS sample. While 
[Anderson et al.l (|2012l ) i s based on t he co nstraints inferred 
from the BAO signal, iReid et al.1 (|2012l ) deals with the 
analysis of redshift-space distortions. Both of these studies 
present constraints on the quantity 



a = ys {Zm)/ys{Zm), 



(21) 



where j/s'*(2m) is the value corresponding to our fiducial 
cosmology. Dropping the priors on uJh and Ub, we obtain 
the constraint ys{zjn) = 0.0745 ± 0.0014, which implies 
a = 1.015 ±0.019. This result is in good agreeme nt with the 



constraints reported in our companion papers: Rei d et al 
(2012) obtain a = 1.023 ± 0.019, while Anders on et al 
(2012) find a — 1.016 ± 0.017 from the pre-reconstruction 
correlation function, and a post- reconstruction "consensus" 
value between ^(s) and P(fc) of a = 1.033 ± 0.017. This 
agreement is a clear demonstration of the consistency be- 
tween the different analysis techniques implemented in these 
studies. 



5 COSMOLOGICAL IMPLICATIONS 

In this Section, we perform a systematic study of the con- 
straints placed on the values of the cosmological parameters 
described Section [4.21 In Section TS. II we present the results 
for the simple ACDM cosmological model with six free pa- 
rameters. In Section [5.21 we discuss our constraints on non- 
flat models. Section 15.31 deals with the constraints on the 
fraction of massive neutrinos. In Section 15.41 we allow for 
non-zero tensor modes. In Section [5.51 we focus on the con- 
straints on the nature of dark energy. Models where the dark 
energy equation of state is constant over time are analysed 
in Section [5. 5. II while Section [5.5.21 explores the constraints 
on the redshift dependence of wue, parametrized according 
to equation (|12|l . Finally, Section [5.5.31 shows the impact of 
allowing also for models with ^2^ 7^ on the constraints 
on Woe- Tables [^niA7l in Appendix \K\ summarize the con- 
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Figure 5. The marginalized constraints in the fim — h plane for 
the ACDM parameter set. The dashed lines show the 68 and 95 
per cent contours obtained using CMB information alone. The 
solid contours correspond to the results obtained from the com- 
bination of CMB data plus the shape of the CMASS ^(s). 



straints obtained in these parameter spaces using different 
combinations of the datasets described in Sections [2] and |3l 



5.1 The ACDM model 

In this Section we focus on the ACDM model and discuss 
the constraints on the parameter space of equation (jlOp . 
The CMB data alone are able to provide tight constraints 
on this parameter space, especially on quantities such as cjb, 
9 and r, whose constraints show almost no variation when 
other datasets are included in the analysis. However, the 
constraints on other parameters are improved by considering 
additional datasets. 

The dashed lines in Fig. [5] show the two-dimensional 
marginalized constraints in the f2m — h plane obtained using 
CMB data alone. The contours show a degenerac y that fol- 
lows approximately a line of constant Q.mh? (,Percival et al.l 
I2OO2I ). This degeneracy limits the accuracy of the one- 
dimensional constraints on these parameters, which from 
the CMB data alone are fim = 0.266 ± 0.024 and h = 
0.7104z0.020. The solid lines in Fig.[5]show the result of com- 
bining the CMB measurements with the CMASS correlation 
function. The extra information contained in the shape of 
^(s) partially breaks this degeneracy, leading to tighter con- 
straints of fim = 0.282 ± 0.015 and h = 0.696 ± 0.012. The 
dashed line in Fig. [2] corresponds to the best fitting model 
obtained in this case. This model gives an excellent match to 
both the location of the BAO peak and the full shape of the 
CMASS correlation function. On scales s > 80/i~^Mpc, the 
model slightly under-predicts the amplitude of ^(s). Note, 
however, that on these scales the individual points in the 
measurement are correlated. Taking into account the full 
covariance matrix, this model gives x^ = 27 for 32 degrees 



Figure 6. The marginalized posterior distribution in the Qm—^k 
plane for the ACDM parameter set extended to allow for non-fiat 
models. The dashed lines show the 68 and 95 per cent contours 
obtained using CMB information alone. The solid contours corre- 
spond to the results obtained from the combination of CMB data 
plus the shape of the CMASS ^{s). The dotted line corresponds 
to the ACDM model, where Q^ = 0. 



of freedom, providing an excellent fit. This model requires a 
real-space bias factor (i.e., c omputed after accounting for the 
boost factor of lKaiseiJllQTSi ') offer = 1.96± 0.09. This value is 
in excellent agreement with the results of (jNuza et al.ll2012l '). 
who estimated a bias factor of &r — 2 from an abundance 
matching analysis of the small and intermediate scale clus- 
tering of the CMASS sample based on the Multidark simu- 
lation. 

The res ults presented h e re are completely consistent 
with those of lAnderson et al.l ([2013), who explored the cos- 
mological implications of the BAO signal in the CMASS 
correlation function. From the combination of this informa- 
tion with the latest data from the WMAP satellite, they find 
f^m = 0.298 ± 0.017 and h = 0.684 ± 0.013 when the param- 
eter space is restricted to the ACDM model. Although this 
agreement is not surprising, as the two analyses are based 
on the same galaxy sample, it is a clear indication of the 
consistency between the two analysis techniques. 

Although consistent within 1 a, the CMASS correla- 
tion function prefers somewhat higher values of f2m than the 
CMB data. This difference can be traced back to the values 
of ys{zm) obtained from these datasets individually. In the 
ACDM parameter space it is possible to obtain a constraint 
on this quantity on the basis of CMB information alone. 
In this case we obtain ys{zui) = 0.0762 ± 0.0018, while the 
CMASS ^(s) gives ys{zm) = 0.0742 ±0.0014. We wiU return 
to this point in Section |6l where we analyse the clustering 
properties of the NGC and SGC sub-samples separately. 

As can be seen in Figure [5] by preferring higher val- 
ues of f2ni, the CMASS correlation function also leads to 
slightly lower values of the Hubble parameter than in the 
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Figure 7. The marginalized, one-dimensional likelihood distribution of the extensions of the ACDM model explored in Sections 15.21 
through l5.5l Dashed lines indicate the constraints obtained from CMB information only, solid lines correspond to the results of the CMB 
plus the shape of the CMASS ^(s), and the dot-dashed lines show full constraints including also BAO and SN data. 



CMB only case. AlthouRh this val ue is lower than the direct 
measurement of lRiess et al.l l|2011f ). t he difference is not sta - 
tisti cally siRnificant. As discussed in [Anderson et al.l (|2012l ) 
and lMehta et al.l (|2012l ) this difference can be reduced if the 
effective number of relativistic species, A^cff, is allowed to 
deviate from the standard value of A'^ff = 3.04. 

As shown in Table IXT] when the SN and BAO data are 
added to the analysis, the results point towards values of 
f2in similar to those of the CMB-I-CMASS case. Combining 
the information from all these datasets, the recovered values 
of f2m and h are similar to the CMB-I-CMASS results and 
the uncertainties are reduced by 33%. In this case we find 
fim = 0.2846lg;;;g^f and h = 0.6941 ± 0.0081. 

Recent analyses have consistently shown evidence of 
a departure from the scale-invaria nt primordial power 
spectrum of sca l ar fluctuations (Sanchez et al. 120061: 



Spergel et al.ll2007l : lKomatsu et al.i 

l|). Our CMB-I-CMASS constraint on the spectral in- 



0.9620lnnna?, increasing the significance of 



2009 , 2011 ; Keisler et all 
2011 " 

deX IS Us — U.OU^U_Q QQgj^, 

this detection to 4.1 a. This limit is almost unchanged 
when all datasets are considered, in which case we get 
Hs = 0.9613lQQQgQ. The deviation from scale-invariance of 
the primordial power spectrum has important implications, 
as most infiationary mo dels predict that the scalar spectral 
index is less than one IjLindel 120081 ). However, these mod- 
els also predict the presence of non-zero tensor primordial 
fiuctuations. As we will see in Section [5]4l although the con- 
straints on ns become weaker when the tensor-to-scalar ra- 
tio, r, is allowed to vary, we also detect a deviation from 
scale invariance at the 99.7 per cent confidence level (CL) 
in this case. 

The results from our study show that the standard 
ACDM model is able to accurately describe all the datasets 
that we have included in our analysis and that the values of 
its basic parameters are constrained to an accuracy higher 
than 5 per cent. In the following sections we focus on con- 
straining possible deviations from this simple model. 



5.2 Non-flat models 

In this Section we drop the assumption of a flat Universe 
and allow for models where ^k 7^ 0. This parameter space is 



poorly constrained by the CMB data due to the so-called ge- 
ometrical degeneracy IjEfstathiou fc Bondlll999l ) relating the 
physical size of the sound horizon at recombination ra(2:«), 
and the angular diameter distance Da{z*). The former de- 
termines the true physical scale of the acoustic oscillations, 
while the later controls its mapping onto angular scales in 
the sky. Models with the same value of O = rs(2,)/DA(2.) 
predict the same position of the acoustic peaks in the CMB 
spectrum and cannot be distinguished on the basis of the pri- 
mary CMB fluctuations alone. This degeneracy is shown by 
the dashed lines in Fig. [S] which correspond to the 68 and 95 
per cent CL contours in the f2m — ^fc plane obtained from the 
CMB data. The dashed line in panel (a) of Fig. [7] shows the 
corresponding marginalized constraints on Q.k, which allow 
for significant deviations from the ACDM model value. In 



this case we obtain fit 



-0.014l°;°^2 and f7„ 



0.32 



+0.10 



As shown by the solid lines in Fig. [6l the constraints 
on j/s(2m) and A{zm) provided by the CMASS correlation 
function are very effective at breaking this degeneracy, lead- 
ing to a drastic decrease in the range of allowed values 
for these parameters. The solid line in panel (a) of Fig. [7] 
corresponds to the posterior distribution of Q.k obtained 
from the CMB-I-CMASS combination, which is in much 
closer agreement with a flat universe. In this case we ob- 
tain fim = 0.285 ± 0.015 and VLk = -0.0043 ± 0.0049. 



[Anderson et al.l (|2012l ) explored the same parameter 
space using the CMASS BAO signal. From the combi- 
nation of this measurement with WMAP data they find 
fim = 0.299 ± 0.016 and fifc = -0.008 ± 0.005. These con- 
straints are in good agreement with findings reported here, 
although they show a preference for slightly higher values of 
the matter density parameter. 

The inclusion of the SN and BAO datasets does not 
significantly improve the results over those obtained using 
the CMB-I-CMASS combination, with a final constraint of 
Qk ~ —0.0045 ± 0.0042 obtained from the combination of 
all datasets. This means that current observations restrict 
possible variations in the spatial curvature of the Universe 
up to a level of An*, ~ 4 x 10"^. 
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5.3 Massive neutrinos 

In the standard ACDM scenario the dark matter compo- 
nent is given entirely by cold dark matter. However, over the 
last decade a number of experiments have shown clear evi- 
dence of neutrino oscillations, implying that the three known 
types of neutrino have a non-zero mass and contribute to 
the total energy budget of the Universe. These observations 
are only sensitive to the mass-squared differences between 
neutrino flavours rather than on their absolute masses. Ab- 
solute neutrino mass measurements can be obtained from 
tritium /3-decay experiments, which at present p rovide up- 
per l i mits of y^ nil, < 6 eV at the 95 per cent CL (jLobashevI 
l2003l : lEitel2005l : lLesgourgues fc Pastoij|2006l ). Future exper- 
iments like KATRIN are e xpected to improve these b ounds 
by an order of magnitude (jOtten &: Weinheimenl2008l ). Un- 
til then, the best observational window into neutrino masses 
is provided by cosmological observations, i n particular by 



the comb ination of CMB an d LSS dat asets JHu et al.l 
Elgar0v et al. 200J; iHannestad 2002; ISanchez et al 



1998; 



20061 : 



Reid et al.ll2010bl : IdePutter et al.ll2oT3 ). A variation in the 



neutrino mass can alter the redshift of matter-radiation 
equality, thereby affecting the CMB power spectrum. Ad- 
ditionally, until the time when they become non-relativistic, 
neutrinos free-stream out of density perturbations, suppress- 
ing the growth of structures on scales smaller than the hori- 
zon at that time, which is a function of their mass. This 
affects the shape of the matter power spectrum and the cor- 
relation function. 

In this Section we explore the constraints on the neu- 
trino fraction, fv. As current estimates of the differences in 
the neutrino mass hierarchy are an order of magnitude lower 
than the constraints on ^ rrii, from cosmological observa- 
tions, these are not yet sensitive to the masses of individ- 
ual neutrino eigenvalues; we therefore assume three neutrino 
species of equal mass. The dashed line in panel (b) of Fig. [7] 
corresponds to the constraints on the neutrino fraction ob- 
tained from CMB data alone. In this case, we find fv < 0.11 
at 95 per cent CL. The solid line in the same panel shows the 
effect of including also the information from the shape of the 
CMASS correlation function, which drastically reduces this 
limit to fv < 0.055 at 95 per cent CL. Our results can be 
converted into constraints on the sum of the three neutrino 
masses using equation (|14|) to obtain ^ m^ < 1.4 eV (95 per 
cent CL) in the CMB only case, and Yl "ii/ < 0.61 eV (95 
per cent CL) for CMB data plus the CMASS ^(s). 

Fig. [8] shows the 68% and 95% constraints in the fim — 
fi, plane. As shown by the dashed lines, when CMB data 
alone is considered, allowing for /i/ 7^ leads to significantly 
weaker constraints on Sim with respect to the ACDM case, 
with its range of allowed values increasing by more than a 
factor of two. The information in the shape of the CMASS 
correlation function improves these constraints, leading to 
flm ~ 0.298 ± 0.019, with a similar accuracy to that of the 
/^ = case. 

In a recent analvsis. lde Putter et al.l (|2012l ) explored the 
constraints on ^ nii, from the angular power spectrum of a 
galaxy sample drawn from BOSS- DR8 following t he CMASS 
selection criteria, as measured bv lHo et al.l (|2012l ). From the 
comb ination of this measur ement with WMAP7 informa- 
tion, Ida Putter et al.l (|2012i ') obtained a limit of Yl 'T'-v < 
0.56 eV at 95% CL, which is relaxed to 5^m„ < 0.90 eV 
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Figure 8. The marginalized posterior distribution in the /^ — Cm 
plane for the ACDM parameter set extended by allowing for a 
non-negligible fraction of massive neutrinos. The dashed lines 
show the 68 and 95 per cent contours obtained using CMB in- 
formation alone. The solid contours correspond to the results ob- 
tained from the combination of CMB data plus the shape of the 
CMASS C(s). 



(95% CL) when a more conservative galaxy bias model 
is implemented. The similarity between these limits and 
our CMB-(-CMASS constraint illustrates the power of us- 
ing the full three dimensional clustering information, which 
can compensate for the much larger volum e probed by the 
sample analysed bv lde Putter et al.l l|2012r ). 

Although not directly sensitive to /^, the additional 
information from SN or BAG measurements improves the 
limits on the neutrino fraction by constraining parame- 
ters which are degenerate with this quantity. Combining 
all datasets we obtain /^ < 0.049 and y^ nii, < 0.51 eV a t 
95 per cent CL. In the analysis of Ide Putter et al.l l)2012f ). 
the inclusion of the SN and Ho measurements provided a 
tighter constraint, with "^rrii, < 0.26 eV at 95% CL and 
X] ^■i' < 0.36 eV (95% CL) for the two galaxy bias models 
they analysed. 

An extension of the current analysis to include informa- 
tion from ^(s) on smaller scales, where it is more sensitive to 
the effect of neutrino free-streaming, could help to improve 
the constraints on the neutrino fraction even furthe r . How - 
ever, as pointed out bv lSwanson. Percival fc Lahavl (|2010|), 
effects related to non-linearities and galaxy bias on these 
scales might impose a limitation on the robustness of clus- 
tering measurements as a means to obtain bounds on the 
neutrino mass. For this reason, the constraints on '^ m^ pre- 
sented here should be regarded as conservative, while the full 
constraining power of the CMASS sample on this quantity 
will be explored in future studies. 
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Figure 9. The marginalized posterior distribution in the rts — 
r plane for the ACDM parameter set extended by allowing for 
non-zero primordial tensor modes. The dashed lines show the 68 
and 95 per cent contours obtained using CMB information alone. 
The solid contours correspond to the results obtained from the 
combination of CMB data plus the shape of the CMASS i{s) . The 
dotted line corresponds to the scale-invariant scalar primordial 
power spectrum, with rig = 1. 



Figure 10. The marginalized posterior distribution in the ei — 
€2 plane for the ACDM parameter set extended by allowing for 
non-zero primordial tensor modes. The dashed lines show the 68 
and 95 per cent contours obtained using CMB information alone. 
The solid contours correspond to the results obtained from the 
combination of CMB data plus the shape of the CMASS §(«)• 
The dot-dashed lines correspond to chaotic inflationary models 
with p = 1, 2 and 4, as indicated by the labels. The dotted line 
corresponds to a constant value of A'^ = 60. 



5.4 Tensor modes 

We now extend the parameter space of equation (fTUl) to in- 
clude the tensor-to-scalar amplitude ratio r. This is the pa- 
rameter space most relevant for the study of inflation as the 
most simple inflationary models predic t non-zero p rimordial 
tensor modes (i.e. gravitational waves. iLindeli2008l ). 

Panel (c) of Fig. [7] shows the marginalized constraints 
on r for the cases of CMB data only (dashed lines) and CMB 
plus the CMASS ^(s) (solid lines). The constraints imposed 
on r by CMB information alone are r < 0.21 (95 per cent 
CL). The CMASS correlation function tightens this limit to 
r < 0.16 at the 95 per cent CL. This result is only marginally 
improved by the additional information of the SN and BAO 
datasets to our final constraint of r < 0.15 (95 per cent CL). 
These results show good agreement with the constraints of 
iKeisler et al] l|201ll ), who found r < 0.17 (95 per cent CL) 
from the combination of the same CMB datasets with BAO 
and Ho measurements. 

Fig. [9] shows the likelihood contours in the rig — r plane 
obtained by means of CMB data alone (dashed lines), and 
its combination with the CMASS ^(s) (solid lines). Tensor 
modes contribute to the CMB temperature power spectrum 
only on large angular scales {£ < 400). An increase in the 
value of r can be compensated for by reducing the ampli- 
tude of the scalar modes, thereby maintaining the total am- 
plitude of the temperature fluctuations at a constant level. 
The consequent decrease of power on smaller angular scales 
can be compensated for by increasing in the scalar s pectra l 
index, rig. Although, as discussed in iKeisler et al.l ((20111), 
the information from the small angular scales of the CMB 



fluctuations provided by SPT does a good job at breaking 
this degeneracy, a residual relation between these parame- 
ters limits the accuracy of their marginalized constraints. By 
also including the information from the shape of the CMASS 
correlation function, it is possible to restrict the range of 
allowed values for these parameters even further. In partic- 
ular, this combination allows us to detect a deviation from 
the scale-invariant primordial power spectrum (indicated by 
the vertical dotted line) with ris < 1 at the 99.7% CL, even 
in the presence of tensor modes. This detection has strong 
implications for the inflationary paradigm. 

We can explore the implications of our results in terms 
of constraints on inflati onary models b y analy sing the hori- 
zon flow parameters of ISchwarz et al.l (|200lh . These are a 
hierarchy of parameters describing the evolution of the Hub- 
ble factor during inflation. The first parameter is given by 
ei = —d\nH{N)/dN, where A^ is the number of e-foldings 
before the end of infiation at which our pivot scale crosses 
the Hubble radius during infiation, and the remaining ones 
are defined through the relation 



din I 



tj+i 



dN 



j^l. 



(22) 



The weak energy condition implies that ei > 0, while a 
necessary condition for infiation is ei < 1 (which implies 
a > 0). The slow- roll approximation can be expressed as 
\ej\ <^ 1, for all j > 0. In this limit, these parameters satisfy 
the relations 



r = 16ei, 



(23) 
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= 1 - 2ei - £2. 



(24) 



These relations can be used to translate our constraints on 
Us and r into the ei — £2 plane. Fig. 1 101 shows the constraints 
obtained in this way. Marginalizing over e2, the combina- 
tion of CMB data plus the CMASS correlation function im- 
plies the limit ei < 0.0097 at the 95 per cent CL. These 
datasets strongly favour models with positive values of £2, 
in which inflation will en d naturally with a viol ation of the 
slow- roll approximation (jLeach fc Liddld 120031 ). From the 
CMB-I-CMASS combination we obtain the limit £2 > at 
the 95.8 per cent CL, which is only marginally improved to 
the 97 per cent CL with the inclusion of the SN and addi- 
tional BAO measurements. 

The horizon flow parameters are related to the infla- 
ton potential V and its derivatives with respect to the in- 
flaton field 0. Then, they can be used to constrain which 
typ e of potentials are compatible with the obs e rvations (see 



Liddle 



entials are compatible with tne obs e rvations (see 
fc LeachllioOM iKinnev et aLll2008l : iFinelli et all 



I2OIOI ). As an example, we explore the constraints on a par- 
ticular class of models, that of the chaotic (or monomial) 
inflation, in which the inflationary phase is driven by a 
potential of the form V{(j)) oc 0^. These models predict a 
simple relation between the horizon flow parameters, the 
power-law index, p, an d the number of e-folds, A^, given by 
l|Leach fc Liddlell2003l l: 



£2 = 



N = 



4 

-ei, 

P 



1 



-1 



(25) 
(26) 



The dot-dashed lines in Fig. [in] correspond to chaotic in- 
flationary models with p = 1, 2 and 4, as indicated by the 
labels. As can be seen from equations (|25|) and (|26|) . a given 
value of A'' corresponds approximately to a constant value 
of £2. For the pivot scale considered here, a plausible upper 
limit for the number of e-fol ds is A^ < 60 (jDodelson fc Huil 
l2003l : iLiddle fc LeachI [2003bl ) . corresponding to £2 > 0.017 
(indicated by the dotted line in Fig. llOf) . If we restrict our 
analysis to this region of the parameter space, we see that 
models with p > 2 are strongly disfavoured by the data. In 
fact, the marginalized distribution for p obtained from the 
CMB-I-CMASS combination after applying this prior implies 
a limit of p < 1.2 at the 95 per cent CL, imposing a strong 
constraint on the viable chaotic inflationary models. 



5.5 The dark energy equation of state 

Until now we have assumed that the dark energy component 
corresponds to a cosmological constant, with a fixed equa- 
tion of state specified by wde = — 1. In this Section, we allow 
for more general dark energy models. In Section 15.5. ll we ex- 
plore the constraints on the value of mde (assumed redshift- 
independent). In Section I5.5.2l we obtain constraints on the 
time evolution of this parameter, parametrized according to 
equation (|12p . Section 15.5.31 deals with the effect of the as- 
sumption of a fiat universe on the constraints on uide- 

In these tests we consider mo dels with wqe. < — 1, c orre- 
sponding to phantom energy fsee lCopeland et al.ll2006l . and 
references therein). When exploring constraints on dynam- 
ical dark energy models, these are allowed to cross the so- 
called phantom divide, «;de = —1. In the framework of gen- 




Figure 11. The marginalized posterior distribution in the 
Sim — WIDE plane for the ACDM parameter set extended by in- 
cluding the redshift-independent value of «!de as an additional 
parameter. The dashed lines show the 68 and 95 per cent con- 
tours obtained using CMB information alone. The solid contours 
correspond to the results obtained from the combination of CMB 
data plus the shape of the CMASS 5(s). The dot-dashed lines 
indicate the results obtained from the full dataset combination 
(CMB-|-CMASS-fSN-(-BAO). The dotted line corresponds to the 
ACDM model, where uiqe = ~1- 



eral relativity, a single fiuid, or a single scalar field without 
higher derivatives, cannot cross t his threshold since it would 
become gravitationally unstab le (jFeng et al.l [2005J; IVikmanI 
I2OO5I : lHi]||2005l : IXia et allbOOSJ ). requiring at least one extra 
degree of freedom. However, models with more degrees of 
freedom are difficult to implement in general dark energy 
studies. Here we fo llow the pa r ametr ized post-Friedmann 
(PPF) approach of iFang et al.l (|2008r ). as implemented in 
CAMB, which provides a simple solution to these problems 
for models in which the dark energy component is smooth 
compared t o the dark matter. Alternatively, as proposed by 
I Zhao et al.l ((20051), it is possible to consider the dark energy 
perturbations using a two-field model, with one of the fields 
being quintessence-like and the other one phant om- like (e.g. 
the quintom model proposed in lFeng et al.ll2005l ) without in- 
troducing new internal degrees of freedom. Both approaches 
give consistent results. 



5. 5. 1 Time-independent dark energy equation of state 

In this Section we explore the constraints on the parame- 
ter set of equation (|10|l extended by including the redshift- 
independent value of wuE as an additional parameter. The 
dashed lines in Fig. [TT] show the two-dimensional marginal- 
ized constraints in the Q.m — wde plane obtained from CMB 
data alone. There is a strong degeneracy between these pa- 
rameters along which different models predict the same an- 
gular position for the peaks in the CMB power spectrum. 
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This is analogous to the geometrical degeneracy described 
in Section 15.21 corresponding to models with constant val- 
ues of Q>. This degeneracy leads to poor one-dimensional 



constraints of ujde 



-l.lStoJg and Q.u 



0.248 



+0.093 



The solid lines in Fig.[TT]show the effect of including the 
CMASS correlation function in the analysis. The constant- 
O degeneracy can be partially broken by providing an ad- 
ditional distance constraint. The constraint on 3/s(zm) pro- 
vided by ^(s) breaks the degeneracy between Q.^ and wde, 
tightening the constraints on the dark energy equation of 
state. The impact of including the CMASS correlation func- 
tion on the marginalized constraints on wde can be seen 
in panel (d) of Fig. [7] where the dashed lines correspond to 
the result of the CMB only case and the solid lines the one 
of the CMB+CMASS combination. In this case we obtain 



fin 



0.295 



+0.041 



and wde 



-0.95lg'2j, in good agreement 



with a cosmological constant. 

From the combination of the BAO signal inferred 
from the CMASS P (k) and £,{s) with WMAP data, 
[Anderson et al.l (|2012h obtained the constraints Jim = 
0.323 ± 0.043 and «;de = -0.87 ± 0.24, in good agreement 
with our findings. As in the previous parameter spaces, this 
is a clear indication of the consistency between the two anal- 
ysis techniques. The extra information in the shape of ^(s) 
improves the constraints on the dark energy equation of 
state by ~20% with respect to the BAO only result, in- 
dicating that, at this redshift, most of the information on 
this parameter is obtaine d through the me a surem ent of j/s. 

In a recent analysis, iMontesano et al.l ([2013) used the 
full shape of the power spectrum of a sample of LRGs from 
the final SDSS-II, combined with a compilation of CMB ex- 
periments, t o obta in the constraint wde = —1.02 ± 0.13. 
iMehta et al.l (|2012l) combined the BAO dist ance measure- 
ment derived by iPadmanabhan et al.l l|2012l ) and IXu et al.l 
((2013) from the same galaxy sample with WMAP data, to 
obtain wde = —0.92 ± 0.13. As these measurements are 
based on observations are lower redshifts, which are more 
sensitive to variations in «;de, they provide slightly tighter 
constraints on this parameter than the CMB-I-CMASS com- 
bination. 

Including also the additional BAO data in the analy- 
sis gives similar results to the CMB+CMASS case, with a 
constraint on the dark energy equation of state of wde = 
— 0.9lloii- When the SN data is considered in the anal- 
ysis instead of the BAO, the resulting constraints are in 
better agreement with the standard ACDM value, with 
wde ~ — 1.054J;Qg7g. It is interesting to note that this re- 
sult is mostly driven by the CMASS+SN combination. In 
fact, the combined information from these two datasets pro- 
vides the constraint wwe = —1.04 ± 0.11, independently 
of any CMB data. Our final constraints obtained from the 
combination of all datasets are shown by the dot-dashed 
lines in Fig. Illl corresponding to f^m = 0.281 ± 0.012 and 
«;de = —1.033^0073. This result is in excellent agreement 
with the standard ACDM model value of wde = — 1, indi- 
cated by a dotted line in FigfTT] 



5.5.2 The time evolution of woe 

In the ACDM model, the equation of state parameter is 
characterized by the fixed value iude = —1 at all times. 
A detection of a deviation from this prediction would be a 
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Figure 12. The marginalized posterior distribution in the 
Wo — ^a plane for the ACDM parameter set extended by allow- 
ing for variations on «;de(i)i parametrized as in equation I I12II . 
The dashed lines show the 68 and 95 per cent contours ob- 
tained using CMB information alone. The solid contours cor- 
respond to the results obtained from the combination of CMB 
data plus the shape of the CMASS £,{s). The dot-dashed lines 
indicate the results obtained from the full dataset combination 
(CMB+CMASS+SN-I-BAO). The dotted lines correspond to the 
canonical values in the ACDM model. 



clear signature of the need of alternative dark energy mod- 
els. In this Section, we explore the constraints on the red- 
shift dependence of uide which we parametrize according to 
equation (|12p . 

The dashed lines in Fig. [12] show the two-dimensional 
marginalized constraints in the wo — Wa plane obtained from 
the CMB data alone. This case provides only weak con- 
straints on these parameters, allowing for models where the 
value of wde can vary significantly over time. The inclu- 
sion of the CMASS correlation function reduces this allowed 
region to a linear degeneracy between wo and Wa which 
can still accommodate large deviations from the ACDM val- 
ues, indicated by the dotted lines. At least a third dataset 
is required to obtain more restrictive constraints. In the 
CMB-FCMASS-hSN case, we obtain wo = -1.09 ± 0.11 and 
Wa = 0.1210.4?, that change to wo = -0.95 ± 0.27 and 
Wa = 0.05^0 61 i^ the SN data are replaced by the addi- 
tional BAO measurements. The dot-dashed lines in Fig. [12] 
correspond to our tightest constraints, obtained by combin- 
ing all datasets, where we obtain the marginalized values 
wo = -1.08 ± 0.11 and Wa = 0.23 ± 0.42. 

A useful quantity to characterize the constraints on the 
redshift evolution of the dark energy equation of state is 
the pivot redshift, Zp, defined as the point where the un- 
certainty on WDE(a) is minimized (JHuterer fc Turnen [200 1[ : 
[Hu fc Jain"2004[: [Albrecht eraL|[200 ^. The parametrization 
of equation H12p implies that this redshift corresponds to the 
scale factor 
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Figure 13. The marginalized posterior distribution in the «>de — 
Hfe plane for the ACDM parameter set extended by allowing for 
simultaneous variations on w)de (assumed time-independent) and 
f2j.. The dashed lines show the 68 and 95 per cent contours ob- 
tained using CMB information alone. The solid contours cor- 
respond to the results obtained from the combination of CMB 
data plus the shape of the CMASS S,{s). The dot-dashed lines 
indicate the results obtained from the full dataset combination 
(CMB+CMASS-fSN-l-BAO). The dotted lines correspond to the 
values of these parameters in the ACDM model. 



Op = 1 -f 



{5wl) 



(27) 



The corresponding pivot redshift for the CMB-I-CMASS 
combination is given by z^ = 1.21, for which we obtain 
Wde(Zp = 1-21) = -0.94 ± 0.20. The pivot redshift for the 
combination of all datasets is Zp = 0.23, which corresponds 
to our tightest constraint on the dark energy equation of 
state, with wde(zp = 0.23) = — 1.040±0.072, in good agree- 
ment with a cosmological constant. 



5.5.3 Dark energy and curvature 

We now explore the constraints on the dark energy equa- 
tion of state (assumed time-independent) when the assump- 
tion of a flat Universe is dropped. This parameter space 
presents similar characteristics to the one studied in Sec- 
tion l5.5.2] where the dark energy equation of stat e is all owed 
evolve over time. As discussed bv lKomatsu et all (|2009l ') and 
I Sanchez et al.l ((20091), when both wde and Qk are allowed 
to vary, the one-dimensional degeneracies corresponding to 
constant values of obtained from the CMB observations 
in the analyses of Sections 15. 21 and lS. 5. II gain an extra degree 
of freedom. As shown by the dashed lines in Fig. 1131 when 
projected in the wue — fife plane, this two-dimensional de- 
generacy extends over a large region of the parameter space. 
The solid lines in Fig. ll3l show the resulting constraints from 
the CMB-fCMASS combination. Although the constraint 
on j/s(2in) provided by the CMASS correlation function sub- 



Figure 14. Large-scale correlation function of the NGC (circles) 
and see (squares) CMASS sub-samples. The dashed line corre- 
sponds to the best-fitting ACDM model obtained by combining 
the CMB data with the information from the shape of the NGC 
correlation function. Although the two measurements exhibit the 
same broad-band shape, in the SGC data the BAO peak has a 
larger amplitude and is located at smaller scales than in the NGC 



stantially reduces the allowed region for these parameters, 
the remaining degeneracy between them corresponds to poor 
one-dimensional marginalized restrictions. 

The distance measurements provided by the additional 
BAO or SN datasets can break the remaining degeneracy, 
leading to meaningful constraints on these parameters. The 
dot-dashed lines in Fig. [13] correspond to the constraints 
obtained with the combination of all four datasets, showing 
good agreement with the ACDM model values (indicated by 
the dotted lines). In this case, we obtain flk ~ —0.0054 ± 
0.0044 and wue = —1.060 ± 0.075, with similar accuracies 
to the constraints obtained when each of these parameters 
are varied independently (Sections 15.21 and I5.5.ip . 



6 THE CLUSTERING SIGNAL IN THE 

NORTHERN AND SOUTHERN GALACTIC 
HEMISPHERES 

Our analysis is based on the full CMASS sample, combin- 
ing the NGC and SGC data. Compared to the NGC, the 
SGC observations correspond to a region with larger av- 
erage Galactic extinction and were taken under higher air- 
mass and sky background and over different periods of time. 
These differences make the NGC-SGC split a sensible cut to 
study the clustering properties of these sub-samples individ- 
ually. In fact, when analysed separately, the clustering of the 
NGC and SGC CMASS sub-samples presents some intrigu- 
ing differences. This can be seen in Fig. 1141 which shows the 
measurements of ^(s) in these two regions, obtained as de- 
scribed in Section (2] It is clear that, although they exhibit 
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the same overall shape, the BAO feature in the SGC has 
a higher amplitude, and its centroid is located at smaller 
scales than in the NGC. In this Section, we explore the sig- 
nificance of these differences and their implications on the 
obtained cosmological constraints. 

IRoss et al.1 (|2012r ) performed a comprehensive analysis 
of the differences between the NGC and SGC CMASS sub- 
samples an d found no treatment of t he data that could alle- 
viate them. ISchlaflv et~aLl (I2OI0I ') and lSchlaflv fc Finkbeineij 
l|2011f ) found small systematic offsets between the colours of 
SDSS objects in the NGC and SGC, which lead to slightly 
different selection crit eria for the CMAS S sample in the two 
galactic hemispheres. I Ross et al.l l|201ll 'l found a 3.2% dif- 
ference in the number density of CMASS targets between 
th e NGC and SGC, whic h redu ces to 0.3% when the offset 
of ISchlaflv fc Finkbeineij (|2011f ) is applied to the galaxies 
in the S GC before applyin g the CMASS selection criteria. 
However, IRoss et all l|2012l ) found that these factors do not 
produce a measurable effect on the clustering signal of the 
SGC CMASS sample, and the differences between the cor- 
relation function of the SGC and NGC remain the same. 

The consistency between the measurements in the NGC 
and SGC can be assessed by examining the difference 
Cngc(s) — ^SGc(s). As these regions correspond to well 
separated volumes, we can neglect the covariance between 
them and estimate the covariance matrix for this differ- 
ence simply as Cdifi = Cngc + Cscc, where Cngc and 
CsGC correspond to the covariance matrices of the individ- 
ual NGC and SGC regions. The consistency of the difference 
^NGc(s) — Cngc(s) with cosmic variance can be estimated 
from its x^ value, with respect to Cdiff. In the range of 
scales used in our analysis, 40 < s/(/i~^Mpc) < 200, we 
find x^ = 53.9 for 41 data points. This number changes to 
X^ = 25.2 for 15 points if the test is restricted to the range 
of scales of the BAO peak (70 < s/(/t~^ Mpc) < 130) . Using 
a different bin size of As = 7/i~^Mpc, IRoss et al.1 (|2012l ') 
performed the same test and found similar values of x'^ per 
degree of freedom. This result shows quantitatively that the 
general shapes of these measurements are in agreement, and 
the differences between them are localized at the scales of 
the acoustic peak. Note, however, that this is the range of 
scales from where the constraints on j/s are obtained. 

Another view of this is presented in Fig. 1151 which 
shows the two-dimensional marginalized constraints in the 
A(zni) — ys(zm) plane. While the two measurements point to- 
wards consistent values oi A^z-m), with A{zm) = 0.426±0.021 
and A{z^) = 0.447 ±0.030 from the NGC and SGC, respec- 
tively, the different locations of the acoustic peak inferred 
from these regions lead to ys(2m) = 0.0762 ± 0.0015 and 
J/s(zm) ~ 0.0704±0.0017, which are approximately 2 a apart. 
Despite the fact that the errors in the SCG correlation func- 
tion are almost a factor two larger than those of its NGC 
counterpart, the accuracies of the constraints on j/s obtained 
from these measurements are similar. This is due to the high 
amplitude of the BAO bump in the SGC £,{s) which, as can 
be seen in Fig. 1141 gives a precise determination of the cen- 
troid of the peak, leading to a slightly smaller than expected 
uncertainty on j/s. As was pointed out in Section [S7T1 within 
the ACDM parameter space the CMB data alone is suffi- 
cient to obtain the estimate ys(zm) = 0.0762 ± 0.0018. This 
value shows a remarkable consistency with the result ob- 
tained from the NGC. A comparison of Figs. [2] and [14] shows 
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Figure 15. The marginalized posterior distribution in the 
A{zYn) —ysizm) plane obtained from the correlation function of the 
NGC (solid lines) and SGC (dashed lines) CMASS sub-samples. 
The contours correspond to the 68 and 95 per cent CL. While 
the two measurements point towards consistent values of A{zin), 
their preferred values of yB{zm) deviate by approximately 2 a. 



that, although the correlation function of the full CMASS 
sample is dominated by that of the NGC, which covers a 
larger volume, adding the data from the SGC moves the 
BAO peak towards somewhat smaller scales, leading to the 
resuh y^izm) = 0.0742 ± 0.0014. 

The conclusion from the tests of IRoss et al.l ([2013) is 
that the differences between the NGC and SGC are simply 
due to a statistical fluctuation. However, as the data in the 
NGC covers a volume 3.7 times larger, providing a better 
knowledge of the survey selection function, for complete- 
ness we also discuss here the constraints on the parameter 
spaces of Section 14.21 obtained from the combination of the 
correlation function of the NGC sub-sample with our CMB 
dataset. We do not consider here, however, the extension 
of the ACDM parameter space in which both wue and i}^ 
are allowed to float since, as discussed in Section [5.5.31 the 
combination of CMB data with a measurement of ^(s) is 
not enough to break the strong degeneracy between these 
parameters. The complete lists of parameter constraints ob- 
tained from the CMB+NGC combination is summarized in 
Table El of Appendix S 

For the ACDM parameter space, the mean values for 
the cosmological parameters obtained in the CMB-I-NGC 
case are in closer agreement with those obtained by means 
of the CMB data alone than in the full CMB-I-CMASS case. 
For example, in the CMB-I-NGC case we flnd constraints 
of n,n = 0.265 ± 0.014 and h = 0.711 ± 0.012, in excellent 
agreement with the CMB only resuhs of Qm = 0.266 ±0.024 
and h = 0.710±0.20. The slightly higher value of the Hubble 
parameter obtained in this case reduces the difference with 
the measurement of iRiess et al.l (|201H ) to the 1 a level. 

When the curvature of the Universe is included as a 
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free parameter, the value of ysiz-m) from the NGC breaks 
the geometrical degeneracy in the CMB data closer to the 
locus of the fiat models, yielding a constraint of Q,k = 
-0.0002 ± 0.0049, completely consistent with the flat Uni- 
verse prediction from the inflationary paradigm. 

When constraining the fraction of massive neutrinos, 
the CMB+NGC combination yields /^ < 0.044 and Y^m^ < 
0.52 eV at 95% CL. These hmits are shghtly tighter 
than those obtained in the CMB+CMASS case. Regard- 
ing the constraints on the tensor-to-scalar ratio, from the 
CMB-I-NGC combination we find r < 0.17 at the 95 per 
cent CL, which is equivalent to the limit found using the 
full CMASS £,{s), albeit with a preference for lower matter 
density values, with fim = 0.276 ± 0.016 

The results for the dark energy related parameter spaces 
also change when the full CMASS ^(s) is replaced by the one 
of the NGC. In this case we obtain weaker constraints, with 
fim = 0.246*^^*12 and wde = -1.14 ± 0.26. When Equa- 
tion (|12p is used to explore the redshift dependence of the 



dark energy equation of state, we find wo — — 1.2ltQg'^ and 
Wa ~ O.l^io and a constraint of woEizp) — —1.21 ±0.26 at 
the pivot redshift of Zp = 0.96. 

In all cases analysed, when we restrict our analysis 
to the NGC-CMASS sub-sample the constraints change 
at most by 1 a. T his is in agreement with the results 
of IRoss et al.l (|2012l ). who found the same level of consis- 
tency. In general, we find that the NGC data points towards 
slightly lower values of Jim and higher ones of h than those 
obtained from the full CMASS sample and in closer agree- 
ment with the CMB only case. It should be emphas ised, 
however, that the extensive tests of IRoss et al.l ((2013), to- 
gether with our internal investigations, show no reason for 
preferring the measurements from the NGC alone to the 
measurements from the full CMASS sample, which provides 
our best picture of the clustering of galaxies at 2: ~ 0.57. 



7 CONCLUSIONS 

We have presented an analysis of the cosmological implica- 
tions of the monopole of the redshift-space two-point cor- 
relation function, ^(s), measured from BOSS-DR9 CMASS 
sample. The large volume and average number density of 
this sample make it ideally suited for large-scale structure 
analysis. The information contained in the full shape of the 
CMASS ^(s) allowed us to obtain accurate constraints of 
the parameters j/s(2m) and A(zin), given by equations 1)19^ 
and (|20p . By adopting an explicit, perturbation-theory based 
model for the correlation function in the mildly non-linear 
regime, and marginalizing over its uncertain parameters, we 
are able to exploit information beyond that in the scale of 
the BAG peak alone. We combined this information with 
that of additional cosmological probes, including CMB, SN, 
and BAG measurements from other data sets, to derive con- 
straints on cosmological parameters. We studied the param- 
eters of the ACDM parameter space, and a number of its 
extensions. The main results from our analysis can be sum- 
marized as follows: 

(i) Our results show that the simple ACDM model is able 
to describe all the datasets that we have included in our anal- 
ysis. Given the different nature of these observations and the 



range of redshifts they probe, this is not a minor achieve- 
ment. The basic parameters of this model are constrained 
to an accuracy better than 5%; a clear demonstration of 
the constraining power of observations in the current era of 
precision cosmology. 

(ii) Fig. [7] summarizes our constraints on possible exten- 
sions of the standard ACDM model. We considered non-flat 
models, massive neutrinos, non-zero primordial tensor fluc- 
tuations, and more general dark energy models. In all of 
these cases the inclusion of the CMASS ^(s) in the analysis 
significantly improves the obtained constraints with respect 
to those obtained using the CMB data alone. Our results 
show no significant evidence of deviations from the ACDM 
picture, which can still be considered as our best cosmolog- 
ical model. 

(iii) The information provided by the CMASS correla- 
tion function is essential to obtain tight constraints on 
the curvature of the Universe. We obtain the constraint 
fife = ~0.0043t°°otl from the CMB+CMASS combination 
which is not significantly improved by adding information 
from SN or other BAO data. 

(iv) When massive neutrinos are considered in the anal- 
ysis, we find a constraint of fv < 0.056 at the 95% CL, 
implying a limit of "^ rriv < 0.61 eV on the sum of the three 
neutrino species. This limit is improved to ^ ruv < 0.51eV 
when the SN and BAO data are added to the analysis. 

(v) When considering tensor modes the CMB-I-CMASS 
combination allowed us to obtain a limit on the tensor-to- 
scalar amplitude ratio of r < 0.16 at the 95% CL, which is 
almost unchanged by considering additional datasets. The 
combination of CMB data with the shape of the CMASS 
correlation function reveals a clear signature of a deviation 
from scale-invariance, with Ws < 1 at the 99.7% CL, also in 
the presence of tensor modes. 

(vi) We explored models where the dark energy com- 
ponent does not correspond to a cosmological constant 
and found no signature of a deviation from the standard 
ACDM model. When the value of woe, assumed time- 
independent, is allowed to vary, the CMB-I-CMASS combi- 
nation provides the constraint wde = — 0.95lg'2o. Interest- 
ingly, the CMASS-f SN combination alone provides a tighter 
constraint, with iude ~ —1.04 ± 0.11, independently of any 
CMB data. Our tighter constraints are obtained from the 
combination of all datasets, with wde ~ —1. 033^0074' ^^ 
good agreement with a cosmological constant. This result 
does not change significantly if the assumption of a flat uni- 
verse is relaxed. We also find no evidence of a redshift evo- 
lution of IfDE. 

(vii) Our results ar e in e xcellent agreemen t with those of 
lAnderson et al.l l|2012l ') and lReid etlol l|2012l '). who explored 
the cosmological implications of the BAO and redshift-space 
distortions measurements in the CMASS sample. This high- 
lights the consistency between the different analysis tech- 
niques implemented in each of these studies, and provides a 
reassuring demonstration of the robustness of our results. 

(viii) We studied the clustering of the NGC and SGC re- 
gions separately. The overall shapes of the correlation func- 
tions in these two sub-samples show good agreement, but 
they differ in the location and amplitude of the BAO peak. 
This transl ates into const raints of j/s(zm) which differ at the 
2 a level. iRoss et al.l (12012) performed a detailed analysis of 
the clustering signal in these regions and found no evidence 
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of additional systematic effects in the SGC data, indicating 
that the observed differences are simply due to a statistical 
fluctuation. For completeness, we explored the constraints 
obtained when the NGC correlation function is used in com- 
bination with our CMB datasets. In all cases our results re- 
main unchanged within 1 a, with the NGC data pointing 
towards slightly lower values of r2m and higher ones of h 
than those obtained from the full CMASS sample. 

The current analysis is based on the first spectroscopic 
data release of BOSS. The larger volume that will be probed 
by subsequent data releases, plus the extension of the anal- 
ysis to the lower redshift BOSS galaxies, will reduce the un- 
certainties in the measurement of ^(s) and the calibration 
of the corrections for potential systematic effects, provid- 
ing even more accurate views of the LSS clustering pattern 
in the Universe. This improvement will be accompanied by 
the release of the CMB power spectra measurement from 
the Planck satellite in early 2013. The combination of these 
datasets will undoubtedly provide new, more stringent con- 
straints on cosmological parameters, and open up the possi- 
bility to explore additional extensions to the ACDM model 
which have not yet been fully explored. 



ACKNOWLEDGEMENTS 

We would like to thank Ryan Keisler for his help with the 
implementation of the SPT likelihood code in COSMOMC. 
We also thank Bradford Benson for helping us find a bug 
in our modifications to CoSMOMC by pointing out a dis- 
crepancy between our constraints on the equation of state 
of dark energy using only CMB data and the reported val- 
ues by the WMAP team. AGS would like to thank all users 
of the Pan-STARRS cluster in Garching for their patience 
and support. CGS and JAR-M acknowledge funding from 
project AYA2010-21766-C03-02 of the Spanish Ministry of 
Science and Innovation (MICINN). AJR is grateful to the 
UK Science and Technology Facilities Council for financial 
support through the grant ST/1001204/1. WJP is grateful 
for support from the the UK Science and Technology Facili- 
ties Research Council, the Leverhulme Trust, and the Euro- 
pean Research Council. FP acknowledges support from the 
Spanish MICINNs Consohder grant MultiDark CSD2QQ9- 
00064. JAR-M is a Ramon y Cajal fellow of the Spanish 
Ministry of Science and Innovation (MICINN). 

Funding for SDSS-III has been provided by the Alfred 
P. Sloan Foundation, the Participating Institutions, the Na- 
tional Science Foundation, and the U.S. Department of En- 
ergy. 

SDSS-III is managed by the Astrophysical Research 
Consortium for the Participating Institutions of the SDSS- 
III Collaboration including the University of Arizona, the 
Brazilian Participation Group, Brookhaven National Lab- 
oratory, University of Cambridge, Carnegie Mellon Uni- 
versity, University of Florida, the French Participation 
Group, the German Participation Group, Harvard Univer- 
sity, the Instituto de Astrofisica de Canarias, the Michigan 
State/Notre Dame/ JIN A Participation Group, Johns Hop- 
kins University, Lawrence Berkeley National Laboratory, 
Max Planck Institute for Astrophysics, Max Planck Insti- 
tute for Extraterrestrial Physics, New Mexico State Univer- 
sity, New York University, Ohio State University, Pennsyl- 



vania State University, University of Portsmouth, Princeton 
University, the Spanish Participation Group, University of 
Tokyo, University of Utah, Vanderbilt University, University 
of Virginia, University of Washington, and Yale University. 
We acknowledge the use of the Legacy Archive for Mi- 
crowave Background Data Analysis (LAMBDA). Support 
for LAMBDA is provided by the NASA Ofiice of Space Sci- 



REFERENCES 

Abazajian K., et al. 2009, ApJS, 182, 543 
Adelman-McCarthy J. K., et al., 2008, ApJS, 175, 297 
Aihara H.,et al., 2011, ApJS, 193, 29 
Albrecht A., Bernstein G., Cahn R., et al., 2006, astro. 



arXiv:astro-ph/0609591 



Anderson L., et al., 2012, submitted to MNRAS, 

arXiv: 1203.6594 
Angulo R., Baugh C. M., Frenk C. S., Bower R. G., Jenkins 

A., Morris S. L., 2005, MNRAS, 362, L25 
Angulo R. E., Baugh C. M., Frenk C. S., Lacey C. G., 2008, 

MNRAS, 383, 755 
Beutler F., Blake C, CoUess M., et al., 2011, MNRAS, 416, 

3017 
Blake C, Glazebrook K., 2003, ApJ, 594, 665 
Blake C, Kazin E. A., Beutler F., et al., 2011, MNRAS, 

418, 1707 
Blanton M. R., Lin H., Lupton R. H., Maley F. M., Young 

N., Zehavi I., Loveday J., 2003, AJ, 125, 2276 
Bolton A. S., et al., in prep. 

Cabre A., Gaztanaga E., 2009, MNRAS, 393, 1183 
Chevallier M., Polarski D., 2001, IJMPD, 10, 213 
Cole S., Percival W. J., Peacock J. A., et al., 2005, MNRAS, 

362, 505 
CoUess M. et al. (The 2dFGRS Team), 2001, MNRAS, 328, 

1039 
Colless M. et al. (The 2dFGRS Team), 2003, 



astro-ph/030658l| 



Conley A., Guy J., Sullivan M., et al., 2011, ApJS, 192, 1 
Copeland E. J., Sami M., Tsujikawa S., 2006, Int. J. Mod- 
ern Phys. D, 15, 1753 
Crocce M., Scoccimarro R., 2008, PhRvD, 77, 023533 
Crocce M., Scoccimarro R., 2006, PhRvD, 73, 063519 
de Putter R., Mena O., Giusarma E., et al., 2012, arXiv, 

arXiv:1201.1909 
Dodelson S., Hui L., 2003, Phys. Rev. Lett., 91, 131301 
Drinkwater M. J., et al., 2010, MNRAS, 401, 1429 
Efstathiou G., Bond J. R., 1999, MNRAS, 304, 75 
Eisenstein D. J., Hu W., 1998, ApJ, 496, 605 
Eisenstein D. J., Seo H.-J., Sirko E., Spergel D. N., 2007, 

ApJ, 664, 675 
Eisenstein D. J., Weinberg D. H., Agol E., et al., 2011, AJ, 

142, 72 
Eisenstein D. J., Zehavi I., Hogg D. W., et al., 2005, ApJ, 

633, 560 
Eitel K., 2005, NuPhS, 143, 197 
Elgar0y 0., Lahav O., Percival W. J., et al., 2002, PhRvL, 

89, 061301 
Fang W., Hu W., Lewis A., 2008, PhRvD, 78, 087303 
Feng B., Wang X., Zhang X., 2005, Physics Letters B, 607, 

35 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological implications of the BOSS-CMASS ^(s) 21 



Finelli F., Hamann J., Leach S.M., Lesgourgues J., 2010, 

JCAP, 04, Oil 
Frieman J. A., Turner M. S., Huterer D., 2008, ARA&A, 

46, 1 
Gaztanaga E., Cabre A., Hui L., 2009, MNRAS, 399, 1663 
Gott J. R., Choi, Y-Y., Park C., Kim J., 2009, ApJ, 695, 

45 
Gott J. R., Slepian Z., 2011, MNRAS, 416, 907 
Gunn J. E., et al. 1998, AJ, 116, 3040 
Gunn J. E., et al. 2006, AJ, 131, 2332 

Guo H., Zehavi I., Zheng Z., 2011, arXiv. larXiV:1111.6598l 
Hamilton A. J. S., 1993, ApJ, 417, 19 
Hamilton A. J. S., Tegmark M., 2004, MNRAS, 349, 115 
Hannestad S., 2002, PhRvD, 66, 125011 
Hinshaw G. et al., 2009, ApJS, 180, 225 
Ho S., Cuesta A., Seo H.-J., et al., 2012, arXiv, 

■arXiv:1201.2137 
Hu W., Eisenstein D. J., Tegmark M., 1998, PhRvL, 80, 

5255 
Hu W., Jain B., 2004, PhRvD, 70, 043009 
Hu W., 2005, PRD, 71, 047301 

Huterer D., Turner M. S., 2001, PhRvD, 64, 123527 
Hiitsi G., 2006, A&A, 449, 891 
Jarosik N. et al., 2011, ApJS, 192, 14 
Jones D. H., et al., 2009, MNRAS, 399, 683 
Kaiser, N. 1987, MNRAS, 227, 1 
Kazin E. A., Blanton M. R., Scoccimarro R., et al., 2010, 

ApJ, 710, 1444 
Kazin E. A., Sanchez A. G., Blanton M. R., 2012, MNRAS, 

419, 3223 
Keisler R., Reichardt C. L., Aird K. A., et al., 2011, ApJ, 

743, 28 
Kim J., Park C., Gott J. R., Dubinski J., 2009, ApJ, 701, 

1547. 
Kim J., Park C., Rossi G., Lee S. M., Gott J. R., 2011, 

JKAS, 44, 217 
Kinney W. H., Kolb E. W., Melchiorri A., Riotto A., 2008, 

PRD, 78, 087302 
Komatsu E., Dunkley J., Nolta M. R., et al., 2009, ApJS, 

180, 330 
Komatsu E., 2010, CQGra, 27, 124010 
Komatsu E., Smith K. M., Dunkley J., et al., 2011, ApJS, 

192, 18 
Landy S. D., Szalay A. S., 1993, ApJ, 412, 64 
Larson D., Dunkley J., Hinshaw G., et al., 2011, ApJS, 192, 

16 
Leach S. M., Liddle A. R., 2003, Phys. Rev. D, 68, 123508 
Lesgourgues J., Pastor S., 2006, Phys. Rep., 429, 307 
Lewis A., Challinor A., Lasenby A., 2000, ApJ, 538, 473 
Lewis A., Bridle S., 2002, PhRvD, 66, 103511 
Liddle A. R., Leach S. M., 2003a, MNRAS, 341, 1151 
Liddle A. R., Leach S. M., 2003b, PRD, 68, 103503 
Linde, A. D., 1983, Phys. Lett. B., 129, 177 
Linde A. D., 2008, LNP, 738, 1 
Linder E. V., 2003, PhRvL, 90, 091301 
Lobashev V. M., 2003, NuPhA, 719, 153 
Manera M., et al., 2012, submitted to MNRAS, 

arXiv: 1203:6609 
Maraston C., et al., in prep. 
Masjedi M., Hogg D. W., Cool R. J., et al., 2006, ApJ, 644, 

54 



Masters K. L., Maraston C, Nichol R. C, et al., 2011, 

MNRAS, 418, 1055 
Matsubara T., 2004, ApJ, 615, 573 
Matsubara T., 2008a, PhRvD, 77, 063530 
Matsubara T., 2008b, PhRvD, 78, 083519 
Mehta K. T., Cuesta A. J., Xu X., Eisenstein D. J., Pad- 

manabhan N., 2012, arXiv. arXiv: 1202.0092. 
Meiksin A., White M., Peacock J. A., 1999, MNRAS, 304, 

851 
Montesano F., Sanchez A. G., Phleps S., 2010, MNRAS, 

408, 2397 
Montesano F., Sanchez A. G., Phleps S., 2012, MNRAS in 

press, larXiv: 1107 .4097 
Nuza S. E., Sanchez A. G., Prada F., et al., 2012, arXiv, 

larXiv: 12 02.6057 

Otten E. W., Weinheimer C, 2008, RPPh, 71, 086201 
Padmanabhan N., Schlegel D. J., Seljak U., et al., 2007, 

MNRAS, 378, 852 
Padmanabhan N., Xu X., Eisenstein D. J., Scalzo R., 

Cuesta A. J., Mehta K. T., Kazin E., 2012, arXiv, 

larXiv: 1202.00901 

Peebles P. J., Ratra B., 2003, RvMP, 75, 559 
Percival W. J., Baugh C. M., Bland-Hawthorn J., et al., 

2001, MNRAS, 327, 1297 
Percival W. J., Cole S., Eisenstein D. J., Nichol R. C, 

Peacock J. A., Pope A. C, Szalay A. S., 2007, MNRAS, 

381, 1053 
Percival W. J., Reid B. A., Eisenstein D. J., et al., 2010, 

MNRAS, 401, 2148 
Percival W. J., Sutherland W., Peacock J. A., et al., 2002, 

MNRAS, 337, 1068 
Perlmutter S., et al., 1999, ApJ, 517, 565 
Prada F., Klypin A., Yepes G., Nuza S. E., Gottloeber S., 

2011, arXiv, arXi v:llll. 28891 
Reid B. A., Percival W. J., Eisenstein D. J., et al., 2010a, 

MNRAS, 404, 60 
Reid B. A., Verde L., Jimenez R., Mena O., 2010b, JCAP, 

1,3 
Reid B. A., et al., 2012, submitted to MNRAS, 

larXiv:1203.6641l 

Riess A. G., et al., 1998, ApJ, 116, 1009 
Riess A. G., et al., 2004, ApJ, 607, 665 
Riess A. G., et al., 2009, ApJ, 699, 539 
Riess A. G., et al., 2011, ApJ, 730, 119 
Ross, A. J., et al. 2011, MNRAS, 417, 1350 
Ross A. J., et al., 2012, submitted to MNRAS, 

arXiv: 1203.6499 
Samushia L., Percival W. J., Guzzo L., et al., 2011, MN- 
RAS, 410, 1993 
Samushia L., Percival W. J., RaccaneUi A., 2012, MNRAS, 

420, 2102 
Samushia L., et al., in prep. 
Sanchez A. G., Baugh C. M., Percival W. J., Peacock J. A., 

Padilla N. D., Cole S., Frenk C. S., Norberg P., 2006, 

MNRAS, 366, 189 
Sanchez A. G., Baugh C. M., Angulo R., 2008, MNRAS, 

390, 1470 
Sanchez A. G., Crocce M., Cabre A., Baugh C. M., 

Gaztariaga E., 2009, MNRAS, 400, 1643 
Schlafiy E. F., Finkbeiner D. P., 2011, ApJ, 737, 103 
Schlafiy E. F., Finkbeiner D. P., Schlegel D. J., Juric M., 

Ivezic Z., Gibson R. R., Knapp G. R., Weaver B. A., 2010, 



© 0000 RAS, MNRAS 000, 000-000 



22 A.G. Sanchez et al. 

ApJ, 725, 1175 
Schlegel D., White M., Eisenstein D., 2009, astro, 2010, 314 
Schwarz D. J., Terrero-Escalante C. A., Garcia A. A., 2001, 

Phys. Lett. B, 517, 243 
Scoccimarro R., Sheth R. K., 2002, MNRAS, 329, 629 
Seo H.-J., Eisenstein D. J., 2007, ApJ, 665, 14 
Seo H.-J., Siegel E. R., Eisenstein D. J., White M., 2008, 

ApJ, 686, 13 
Seo H.-J., Eckel J., Eisenstein D. J., et al., 2010, ApJ, 720, 

1650 
Seo H.-J., Ho S., White M., et al., 2012, AAS, 219, #402.02 
Smith R. E., Scoccimarro R., Sheth R. K., 2008, PhRvD, 

77, 043525 
Smee S., et al., in prep. 
Spergel D. N., Verde L., Peiris H. V., et al., 2003, ApJS, 

148, 175 
Spergel D. N., Bean R., Dore O., et al., 2007, ApJS, 170, 

377 
Swanson M. E. C, Tegmark M., Hamihon A. J. S., Hill 

J. C, 2008, MNRAS, 387, 1391 
Swanson M. E. C, Percival W. J., Lahav O., 2010, MN- 
RAS, 409, 3 
Taruya A., Nishimichi T., Saito S., Hiramatsu T., 2009, 

PhRvD, 80, 123503 
Tegmark M., Blanton M. R., Strauss M. A., et al., 2004, 

ApJ, 606, 702 
Tojeiro R., et al, 2012, submitted to MNRAS, 

larXiv:1203.6565l 

Vikman A., 2005, PRD, 71, 023515 

Weinheimer C, 2003, arXiv:he p-ex/0306057| 

White M., et al., 2011, ApJ, 728, 126 

Xia J., Cai Y., Qiu T., Zhao G., Zhang X., (2008), Int. J. 

Modern Phys. D, 17, 1229 
Xu X., Padmanabhan N., Eisenstein D. J., Mehta K. T., 

Cuesta A. J., 2012, MNRAS submitted, arXiv: 1202.0091. 
York D. G., Adelman J., Anderson J. E., Jr., et al., 2000, 

AJ, 120, 1579 
Zehavi I., Blanton M. R., Frieman J. A., et al., 2002, ApJ, 

571, 172 
Zhao G., Xia J., Li M., Feng B., Zhang X., 2005, PRD, 

123515 



APPENDIX A: SUMMARY OF THE 
OBTAINED COSMOLOGICAL CONSTRAINTS 

In this section we summarize the constraints on cosmo- 
logical parameters obtained using different combinations of 
the datasets described in Sections [2] and |3l Table lAll lists 
the 68% confidence limits on the parameters of the ACDM 
model, as discussed in Section [5. II Tables IX21IA7I correspond 
to the extensions of this parameter space analysed in Sec- 
tions 15.21 to 15.51 Finally, Table IA8I presents the constraints 
on these parameter spaces, obtained from the combination 
of the correlation function of the NGC sub-sample with the 
CMB data. 
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Table Al. The marginalized 68% constraints on the cosmological parameters of the ACDM model obtained using different combinations 
of the datasets described in Section 12.11 and [S] 
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Table A2. The marginalized 68% allowed regions on the cosmological parameters of the ACDM model extended by adding Cfe as a 
free parameter, obtained using different combinations of the datasets described in Section |2. II and [3l 
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Table A3. The marginalized 68% allowed regions on the cosmological parameters of the ACDM model extended by adding /j, as a free 
parameter, obtained using different combinations of the datasets described in Section |2. II and l3l 
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Table A4. The marginalized 68% allowed regions on the cosmological parameters of the ACDM model extended by allowing for non-zero 
primordial tensor modes, obtained using different combinations of the datasets described in Section |2. II and l3l 
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Table A5. The marginalized 68% allowed regions on the cosmological parameters of the ACDM model extended by including the redshift- 
independent value of ioqe as an additional parameter, obtained using different combinations of the datasets described in Section |2. II and 

m 





CMB 


CMB + BOSS 


CMB + CMASS 
+SN 


CMB + CMASS 
+BAO 


CMB + CMASS 
+ BAG + SN 


JWde 


1 1C+0.39 
-i-iO_0.39 


q"S+0-2l 
u.yo_Q 20 


1 n^4+o-077 

— i.U04_(, Q7g 


-0.91«:}} 


1 f.qq+O-073 


100 e 


1 0410+00016 


1 0410+00016 


1 0406+00015 


l-0409t°:°°l^g 


1 040^+00015 


100 tJb 


2.219t°:«g 


r, r|r,|-, + 0.043 


r, r,„ +0.039 

^■^11-0.038 


r, 991+0.041 

^■^^1-0.041 


^r,ip,+0.039 

Z.ZlU_g Q3g 


lOOcJdm 


ll-2ltr/7 


11.33«-i? 


ii.48t°;f3 


ii.24t°;f3 


11.58_q;32 


r 


0.0847l°:rv? 


0.083li°:»?g 


0.0819l°:»r4 


0.0840j;°:~?g 


n 0814+0-0057 
'J-'J«^4q QQg3 


Us 


0.9651°;°!} 


0.964t°:°li 


0.9615t°:»^I 


0.966tHii 


9606+00096 


In(lOiOAs) 


Q nci +0.030 


q nH'^+0030 


q nSzl+0-028 
3.U»4_o(,28 


q nsi +0-030 


q 087+0-028 

O.UO(_Q Q2g 


f^DE 


PI 7^.-1+0.088 


0.7041H1? 


0.729lHi^ 


0.702lHrr 


710+0-012 


Hm 


0.248t°;°fg 


0.2951°:°!^ 


0.270l°:°i« 


0.297l°:0rT 


o.28itEJ:S^ 


CT8 


0.86t°;» 


n sm +0.085 


0.842lHi 


0.7871^^4 


0.840tH^^ 


to/Gyr 


13.69l°:^5 


13.82t°:l^ 


13.74l°:°'7^ 


-TQ on+0-085 
l'3-»2_o.084 


13 763+0071 
±o.(UO_Q 072 


■Src 


10.411- 


10.3tl-J 


10.2t};l 


10.3tlJ 


io.2tl:l 


h 


p -7+0.15 
•J- ''-0.14 


n R«'^+0-054 
0.683_Q05g 


(. 710+0.020 

u. ( iO_g Q20 


674+0025 


701+0-016 


£'v(^m)/MpC 


1993t9o 


2045t« 


2018+^^ 


204411° 


2030tlg 


f{Zm) 


n 7P;c;+0.025 


n 7K/1+0.022 


760+0016 


n 748+0-019 

U-'«_Q|3jg 


764+0-015 
U-'04_Q(,j5 



Table A6. The marginalized 68% allowed regions on the cosmological parameters of the ACDM model extended by allowing for variations 
on tooE(a) (parametrized according to equation [T2]l , obtained using different combinations of the datasets described in Section |2. 11 and 

El 





CMB 


CMB + CMASS 


CMB + CMASS 
+SN 


CMB + CMASS 
+BAO 


CMB + CMASS 
+ BAD + SN 


wo 


-, 19+0.52 


-, 1 9+0.61 


-l.09t°;ll 


-0.95«-|? 


1 08+0-11 


Wa 


-o-3tl:? 


0.32l°:°« 


0.121°:!? 


0.05l°:«l 


23+0.42 


1000 


1.0409l°:Ooi6 


1.0409l°:Ooi6 


1 nzina+0-0015 
1.0408_oooig 


1.0409l°:°gl^ 


1.0408l°:~i^ 


lOOwfc 


2.219l°:°g 


2.218l°:°« 


9 91 c:+0-040 
^■^10-0.040 


2.218l°:°°4r 


0.022lt°;°S 


lOOo^dm 


11.22+y^ 


ll-3ltgi^ 


ii.4ot°:« 


11 98+0-48 
11.28_o47 


11 q8+0.47 
11.38_o47 


r 


0.0852l°:»^^ 


0.0833+°oo''6"7'^ 


0.0823l°:°gr7 


O8q'?+0-0061 

U.U«dCS_Q QQgg 


om25+_lZll 


Us 


n Qfic:+0.011 


o.965l°:°ii 


0.963l°:°n 


0.965l°:°li 


0.963l°:°ii 


ln(10iOA,) 


q nsq+0.030 
J.UOJ_Q Q29 


q p,o9+0.030 


q |-,Qq + 0.029 


3.080l°:°29 


q nsq+0.030 


f^DE 


760+0-081 
U.'DU_(j Qg7 


r, 799+0.081 
U.'^^_0.09l 


730+0016 


0.706l°:°g 


724+0.014 


flm 


9qQ+0.087 


n 970+0.091 
'J-^'°-0.081 


269+°°^'' 
u.^oy_Q Qig 


o.294J:°:°l 


976+0-014 


0-8 


0.87t°:ll 


r, 09+0.11 
0-82_o.ii 


n a'i9+0-049 


0.792lg:°r7 


09-1+0.048 


to/Gyr 


13.64+°;i 


13 70+0.16 
^■^■'^-0.16 


13 763+0.089 


13.827t°:°|t 


-iq Qo+0.083 

iJ.OU_Q Qg3 


Zve 


10.4111 


10.3tlJ 


10.2tlJ 


10.3tlJ 


lO.StlJ 


h 


+0.14 
0.78_(, 14 


0.72l°:ll 


7-19+0.020 


n fi8n+0-038 

0-08U_o(,3g 


0.070l°:°i6 


Dv(2™)/Mpc 


1974t|6 


2040j:g 


2027t^^ 


2046^20 


2038±1| 


/(^m) 


f. 7qq+0.077 
U.CdCl_Q Q7g 


0.770i[J:°« 


n -7^0+0.022 
0.7bb_oo22 


7^0+0.040 

U.COCl_Q Q4Q 


+0.019 
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Table A7. The marginalized 68% allowed regions on the cosmological parameters of the ACDM model extended by allowing for 
simultaneous variations on ujoe (assumed time-independent) and fij., obtained using different combinations of the datasets described in 
Section [211 and [3l 





CMB 


CMB + CMASS 


CMB + CMASS 
+SN 


CMB + CMASS 
+BAO 


CMB + CMASS 
+ BAG + SN 


Q.K 


n npfi+ooas 

U.U/0_(,g33 


n nn99+''-oo68 
u.uu^y_g Q(,g4 


nnn^i +0-0048 

-U.UUt>i_Q 0043 


nnni 3+0-0064 


n nn"^4+o-oo44 

U.UUa4_Q (,(,44 


S^DE 


-o-9itg-i? 


-1 n7+''-3'' 

^•'J'-0.38 


-1 n7n+°°™ 


-0.946t°:rg 


1 nfin+O-075 


lOOG 


1 n4in+oooi6 


1 n419+°°016 


1 0411+0-0016 


1 n.11 1+0.0015 


1 041+0-0016 
^•^^4^-0.0016 


fOOcjfc 


r, ri-,0 + 0.041 


2.224lHg 


2.2211HS 


994+0.038 


2.220t°:°S 


fOOcjdm 


11 19+°"'^ 


ii-i6t°:f5 


ll-18t211 


n 11+0-46 


ll-18tg11 


T 


0.0839±Hm 


0.0843tHo66 


0.0844lg:»^i 


0.0850tHolo 


0-OS50l:°oToll 


ris 


0.964lHii 


n ncc+0.011 


n qfi4+o-oio 


0.965lHlg 


0.9641H11 


In(lOiOAs) 


3.078lH^? 


O.U(3_Q Q3Q 


o 079+0-029 


3.08llH^^ 


3.080tH^S 


f^DE 


o.62tg:l^ 


p ~r,c + 0.069 


p -qo+0.017 


0.7071^^7' 


0.7301H11 


nm 


o.4oj:«:- 


277+0-059 


r, riyi +0.016 


294+°°23 
U.ZM4_Q 023 


97^+0.012 


0-8 


0.77l«:i^ 


0.826to;n 


n 0-39+0.036 


0.7951HS 


0.829lHi 


to/Gyr 


i4.7ti:} 


13.90«:- 


10 07+0.23 
icS.y^_(, 23 


13.88+°f 


13.99t°:22S 


■2rc 


lo.sti:^ 


w.stli 


lO.Stl'i" 


10.4+- 


10.4tl;l 


h 


U-O^_0.i6 


n 707+0-087 


n 7n'?+o-02i 


n fi77+0.028 
'J-O"-0.029 


omst°Te 


Dy{Zm) 


2245^2?? 


206lt« 


2054l« 


2053t|? 


206ll|0 


fiZn.) 


n si4+0-0'^3 


n 7fiS+0-041 
U.7b8„oo37 


+0.017 
'J- '"'-0.017 


7^4+0.026 
U.(^t)4_oo26 


n 7fiS+0-015 
'J-'O»-0.015 



Table A8. The marginalized 68% allowed regions on the cosmological parameters of the ACDM model and its extensions, obtained by 
combining the CMB data with the correlation function of the NGC CMASS sub-sample. 





ACDM 


ACDM+Qfe 


ACDM+/^ 


ACDM-l-r 


ACDM-|-«)DE 


ACDM+«)DE(a) 


looe 


l.U411_(, QQj^g 


1 nAI 1+0-0015 
1.0411_(, (,oig 


1 nzii 1+0.0014 

1.0411_oooi4 


1 0411+'^'°015 
l.U4il_Q 0015 


1.0409t°;°»l^ 


1.0409l°;»°l«g 


I^dm 


11 14+0.28 
44-44_o 28 


11 1O+0.46 
tl-t»_0.46 


1 1 QQ'4"0.28 
i_L.zo_Q_28 


n 10+0-28 
ll.lU_o 28 


2.217«:«g 


2.217t°;g« 


100 (.Jb 


o Qon+0.039 

Z.ZZO_Q Q3g 


r, 900 + 0.040 
Z.ZZCl_g Q40 


2.2241H1 


9 904+0.041 
Z.Zd4_o 041 


ll-32tr46 


11.29tO;« 


r 


0.0850«:g°f7 


0-OSi8t'oToll 


0.08621^^76 


n ns/19+0.0060 
0.0842_oo066 


0.0S33+_lZll 


0.0836t°:«r8 


ns 


0.9666t°:«g^^ 


0.966lHn 


n nc7o+0.0095 
0.9b78_ooo95 


971+0-010 

u.y'i_o.oio 


n ncoc+0.0111 
0.9636_o.oili 


0.9641o+«oi°i\" 


In(lOiOAa) 


smot^oll 


o nsn+0-030 


o 077+0.031 
■^■'^"-0.031 


070+0.029 

J.U/y_o 029 


n»9+0.029 
J.USZ_o 029 


3.082t°:gfg 


Hfe 





-0 0002+°°°*3 

U.UUUZ_(, Q(j4g 














u 








< 0.044 (95% CL) 











r 











< 0.17 (95% CL) 








""DE (wo) 


-1 


-1 


-1 


-1 


_1 14+0.26 
J- -^^-0.27 


1 91+0.61 
-1-^1-0.79 


Ula 





-1 


-1 








14+^° 


E^"- 








<0.52eV (95% CL) 











f^DE 


0-735l°:gl* 


7OO+0.015 


790+0.016 
'-'-'^■^-0.017 


0.738l°:°il 


7CO+0.045 
U./OJ_o 042 


0.756tO;«ri 


f^m 


0.265l°:gll 


9^7+0.015 
U-^O' -0.015 


0.2761°:°!^ 


0.262J:H11 


0.246t°;SS 


0.244t0;g?l 


a-8 


n oiq+0.018 


0.814lH^^ 


0.764lHl^ 


n S1 1+0.018 

O.»ii_o.oi8 


0.86ll«:J0 


0.87«:11 


to/Gyr 


13.727«:«?? 


13.74t°;g 


iqsis+o-io 

ici.»i»_o 098 


10 710+0.074 

lo.'lo_o.o75 


13.69j;«:i^ 


nrntrA 


■2^re 


io.4l}:i 


10.4ti;} 


10.5t\i 


10.3ti:? 


io.3tl:i 


w.stli 


h 


711+0.012 


n 7nR+0-0i7 

U- '^'J*<-0.017 


0S9Stlf,l 


r, 710 + 0.013 

U.'1.3_0.013 


7K1 +0.078 
'J- '0^-0.073 


77+0.12 
U.C^_0 12 


I?v(-2m)/MpC 


2005^20 


2010tl^ 


2030l^? 


2OOII22 


19971^^ 


1996^4? 


/(^m) 


n 7t;/i+0.012 
'^•'04-0.012 


n 7/1C+0.013 
U.l^40_ooi3 


(. 7^„+0.014 
U.^OCi_o 014 


740+''-"^2 

U.^4U_y 012 


7^57+0.028 
'J-'S'-0.027 


7f;o+0.064 
U./DU_oo7l 
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